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The Abundances of Light Neutron-Capture Elements in Planetary 
Nebulae — III. The Impact of New Atomic Data on Nebular 
Selenium and Krypton Abundance Determinations 

N. C. Sterling 1 , R. L. Porter 2 , & Harriet L. Dinerstein 3 

ABSTRACT 

The detection of neutron(n)-capture elements in several planetary nebulae 
(PNe) has provided a new means of investigating s-process nucleosynthesis in 
low-mass stars. However, a lack of atomic data has inhibited accurate trans-iron 
element abundance determinations in astrophysical nebulae. Recently, photoion¬ 
ization and recombination data were determined for Se and Kr, the two most 
widely detected n-capture elements in nebular spectra. We have incorporated 
these new data into the photoionization code Cloudy. To test the atomic data, 
numerical models were computed for 15 PNe that exhibit emission lines from 
multiple Kr ions. We found systematic discrepancies between the predicted and 
observed emission lines that are most likely caused by inaccurate photoioniza¬ 
tion and recombination data. These discrepancies were removed by adjusting 
the Kr + -Kr 3+ photoionization cross sections within their cited uncertainties and 
the dielectronic recombination rate coefficients by slightly larger amounts. From 
grids of models spanning the physical conditions encountered in PNe, we derive 
new, broadly applicable ionization correction factor (ICF) formulae for calculat¬ 
ing Se and Kr elemental abundances. The ICFs were applied to our previous 
survey of near-infrared [Kr III] and [Se IV] emission lines in 120 PNe. The re¬ 
vised Se and Kr abundances are 0.1-0.3 dex lower than former estimates, with 
average values of [Se/(0, Ar)] = 0.12±0.27 and [Kr/(0, Ar)] = 0.82±0.29, but 
correlations previously found between their abundances and other nebular and 
stellar properties are unaffected. We also find a tendency for high-velocity PNe 
that can be associated with the Galactic thick disk to exhibit larger s-process 
enrichments than low-velocity PNe belonging to the thin disk population. 
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Subject headings: planetary nebulae: general—nucleosynthesis, abundances— 
stars: AGB and post-AGB—atomic data— infrared: general 


1. INTRODUCTION 


Despite their low cosmic abundances, neutron(n)-capture elements (atomic number Z > 
30) can provide significant insight into the nucleosynthetic histories of stellar populations, 
mixing processes and the internal struct ure of asymptotic giant branch (AGB) stars, and the 


chemical evolution of the Universe (e.g., 

Busso et al 

2001; 

Herwig 

2005; 

Sneden et al. 

2008; 

Karakas et al. 

2009; Simon et al. 2010; 

McWilliam et al. 

2013n. However, much remains 

unknown about the production of n-capture elements (Busso et al. 

1999; 

Sneden et al. 

2008; 

Thielemann et al. 2011; Karakas & Lattanzio 2014: 

Triooclla et al 

20U 

A complicating the 


interpretation of observations of these species. 


Our understanding of n-capture nucleosynthesis has largely been construc t ed from stel- 


lar abundance patterns (A 

derrill 

1952; 

Smith & Lambert 

1990; 

Busso et al. 

1999; 

Travaglio et al. 

2004; 

Sneden et al. 

2008; 

Roederer et al. 

2010, and references therein). In recent years, n- 


capture element emission lines have been detected in a growing number of astrophysical 


nebulae, including 

alanetary nebulae (PNe 

: Liu et al. 

2004b; Zhang et a 

1. 200, 

5; Likkcl et al. 

2006|; 

Sharoee et a 

. 2007; 

Sterling et al. 

20071; Sterling & Dinerstein 

20081 

Sterling et al. 

20Q9|; 

Otsuka et a . 

201(1 

Fang & Liu 201 

1; Otsuka et al. 2011: 

Garcia-Roi 

as et al. 20121: 

Otsuka & Taiitsu 2 

013) 

, H II regions (Asp 

in 19941; Lumsden & Puxlev 

1996; 

Baldwin et al. 

20001; 

Puxlev et al. 

2000 

; Okumura et al. 2( 

101; Blum & McGregor 

2008 

Roman-Lopes et al. 

2009j), and the interstellar medium of other galaxies 

Vanzi et al. 

2008) 

. The large number 


of detections demonstrates the potential of nebular spectroscopy for investigating the abun¬ 
dances and production of n-capture elements. 


Trans-iron element abundances in PNe are of particular interest, since these nuclei can 
be produced via slow n-capture nucleosynthesis (the s-process) in PN progenitor stars. Dur¬ 
ing the thermally-pulsing AGB evolutionary stage, the s-process can occur in the intershell 
region between the H- and He-burning shells. Free neutrons are released predominantly by re¬ 
captures onto 13 C (or 22 Ne in AGB stars with mass > 4-5 M 0 ). Iron-peak nuclei act as seeds 
from which s-process nuclei grow, alternatel y capturing neutrons and undergoing d -decays 
to transform into nuclei of heavier elements ( Busso et al. 19991: Kappeler et al. 2011 b Con¬ 
vective dredge-up events late in the A GB phase transport the s-process enriched material to 
the stellar envelope along with carbon (Busso et al. 1999: Hcrwig 2005: Karakas & Lattanzio 


20140. Stellar winds and PN ejection supply the surrounding interstellar medium with en¬ 


riched material, which is eventually incorporated into new generations of stars. 
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Nebular spectroscopy reveals unique information about s-process nucleosynthesis, pro¬ 
viding access to n-capture elements that are difficult or impossible to detect in AGB stars. 
Elemental abundances in PNe also reveal chemical enrichments after the cessation of nucle¬ 
osynthesis and dredge-up, which is ongoing in thermally-pulsing AGB stars. In addition, 
stars that are challe nging to study during the AGB phase due to dust obscuration or high 
mass- loss rates (e.g., Ventura fe Marigo 20091; Dell’Agli et al. 2015 1 produce readily observed 
PNe ( Sharpee et al. 20071: Sterling V Dinerstein boosi hereafter Paper II). 


Selenium and krypton are the two most widely detected trans-iron elements in PNe 
and other nebulae, owing to their relatively strong near-infrared fine-structur e tra nsitions 
[Kr III] 2.199 and [Se IV] 2.287 //m first identified by Dinerstein ( 2001 1. In Paper II . Se and 
Kr abundances were derived in over 80 PNe from observations of these K band transitions, 
providing the first large-scale investigation of s-process enrichments in PNe. 

However, abundance determinations of Se and Kr — and indeed all n-capture ele¬ 
ments — in astrophysical nebulae are plagued by uncertainties. The most important of 
these uncertainties stems from th e absence of reliable atom i c data . Transition probabil¬ 
ities 


( Biemont fe Hansen 1986allbl. Il987 : 


Biemont et a 


1988 


1995ll and effective collision 


strengths for electron-impact excitation ( Schdningl 1997 : Schoning fe Butler 1998 1 have been 
determined for some n-capture element ions, allowing ionic abundances to be derived. How¬ 
ever, atomic data needed to reliably correct for the abundances of unobserved ionization 
stages are unknown. Such “ionization correction fact ors,” or ICFs, are most reliably derived 
via numerical simulations of astrophysical nebulae (Stasinska 197 si: Kingsburgh fc Barlow 
1994 : Kwitter fe Henry 2001 : Rodriguez fe Rubin 20051: Delgado-Inglada et al. 2014 1. How¬ 
ever, the simulations are only as accurate as the underlying atomic physics that they in¬ 
clude. Specifically, ionization equilibrium solutions in photoionized nebulae rely on accurate 
photoionization (PI) cross sections and rate coefficients for radiative recombination (RR), 
clielectronic recombination (DR), and charge transfer (CT). Until recently, these data were 
unknown for nearly all n-capture element ions. 


In the absence of these data, [Sterling et al . (2007, hereafter Paper I) used approxima¬ 
tions for the necessary atomic parameters to derive IGF prescriptions for Se and Kr. They 
utilized hydrogenic PI cross sections and RR rate coefficients self-consistently derived from 
the PI data, CT rate coefficients of Cu and Zn ions of the same electronic charge, and DR 
rate coefficients equal to the average rate coefficients of C, N, and O ions of the same charge. 
The Kr atomic data was adjusted empirically in order to reproduce emission line strengths 
in PNe in which multiple Kr ions had been detected. Monte Carlo simulations were run to 
estimate the sensitivity of the derived Se and Kr abundances to atomic data uncertainties. 
Those simulations indicated that the poorly-known atomic data can lead to errors in derived 
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Se and Kr abundances by a factor of two or possibly more, if the approximate data are 
substantially far from the actual values. 


Sterling &; Witthoeftl (1201 lh and ISterlingl (120111) recently computed PI cross sections 


and RR and DR rate coefficients for the first six ions of Se and Kr, to address the lack of 
reliable atomic data for thes e elements. The P I cross sections were benchmarked against 


experimental measurements (Lu et al. 


; Eu 


2006t iBizau et al.l l201ll : ISterling et al. 


2011: Esteves et al. 2011: E 

Isteves et al. 

2012h. leading to improved accuracv for some ions 

(Sterling 

2011 

). In addition 

Sterling & Stancil (2011 

) computed CT rate coefficients for low- 


charge ions of several n-capture elements detected in astrophysical nebulae. Uncertainties 
were estimated in all of these calculations, enabling a quantitative analysis of the sensitivity 
of abundance determinations to atomic data uncertainties. We will present the results of 
that investigation in a forthcoming paper. 

To apply the new atomic data to nebular Se and Kr abundance determinations, we have 
added the new PI. RR, DR, and CT data to the widely-used photoionization code Cloudy 
(iFerland et al.l 2013). The purpose of the present study is twofold: we test the new atomic 
data against the observed ionization balance in individual PNe, and derive new ICF formulae 
to account for the abundances of unobserved Se and Kr ions in ionized nebulae. 

To address the first goal, we model individual PNe in which multiple Kr ions have been 
detected and compare the predicted and observed emission line intensities of the detected 
ions. Unfortunately the Se atomic data cannot be tested in a similar manner as only one 
ion, [Se IV], has been unambiguously detected in PNe (with the possible exception of Hb 12; 


Sterl ing et al. 2013). PNe are valuable laboratories for evaluating the accuracy of atomic 


data, owing to their high surface brightness that makes it possible to detect emission lines 


^ -—-o'- - --— ? 

terminations against observed PN spectra (e.g., 

Drake & 

dobbins 19721; 

Garst 

ang et al. 19781; 

Dennefcld & Si 

;asiiiskal 

1983; 

Keenan & Aggarwal 1987: 

Keenan et al. 

1992; 

Bicmont et al. 

1995; 

Liu et al. 

1995; Kastner & Bhatia 1996; Sc 

loning 1997; Kisielius et al. 1998; van Hoof et al. 

2000; 

Chen & Pradhan 

2000; 

Wang et al. 2004; 

Bautista 

2009; McNabb et a 

1. 2013). 


In order to derive new analytical ICFs for Se and Kr, we computed large grids of models, 
spanning the range of physical conditions found in most PNe. Following the methodology of 
Pap er I], we searched for correlations between the fractional ionic abundances of observed Se 
and Kr ions and those of routinely detected light elements. We fitted the correlations with 
functions to derive ICF prescriptions that can be applied to optical and near-infrared obser¬ 
vations of Se and Kr. Using the new ICFs, we derive more accurate Se and Kr abundances 
for the PN survey of Paper II. 
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The structure of this paper is as follows. In §2, we briefly describe our alterations 
of Cloudy to include the new Se and Kr atomic data. In §3, we present detailed models 
of 15 PNe in which multiple Kr ions have been detected, in order to test the new atomic 
data’s ability to reproduce observed line intensities. A large grid of models encompassing 
the range of physical conditions encountered in most PNe is presented in §4, and we derive 
new formulae for Se and Kr ICFs from those results. In §5, w e use the newly-derived ICFs 
to revise the Se and Kr abundance determinations of Paper II in 120 PNe, before providing 
concluding remarks in §6. 


2. THE INCORPORATION OF SELENIUM AND KRYPTON INTO 

CLOUDY 

The modifications described in this section are relative to the C13 release of Cloudy 
(IFerland et al.l 120131) . We emphasize the updated atomic database for Se and Kr, but all 
elements from Ga through Kr were added to the code. Of the elements we added to Cloudy, 
only Se and Kr were activated in the models described in the subsequent sections. 


Atomic weights are taken from the same sourc e as for lighter elemen ts (ICoplenll200ll) . 


The default abundances are the solar values from Asplund et al. 


tion of interstellar medium abundances and dust depletion factors 


(2009), with the excep- 


Savage & Sembach 1996 


Welt v et al . 1999 ). Data for the hydrogenic and He-like isoelectronic sequences were re¬ 
computed to extend to Z = 36. Th ese calculatio n s are s traightforward ap p licatio ns of 
the analytical t r eatme nts described in IPorter et al.l (120051 ). IPorter fc Ferlandl (120071 ). and 
Luridiana et ah (120091) . 


We inco r porat ed the direct valence shell PI cross sections of lSterling &; Witthoeftl (120111) 
and iSterlingj (20111) for the neutral and first five ionization stages of Se and Kr. Following 
the treatment of other species in Cloudy, resonant photoionization was not included due to 
uncertainties in the calculated resonance energies (IVerner et al.1 19961) • For all other Se and 
Kr io ns an d those of Ga, Ge, As, and Br, we adopt hydrogenic PI cross sections as described 


m 


Paper I. 


RR and low-temperatu re D R rate coe f ficien ts for Se°-Se 5+ and Kr°-Kr 5+ are taken from 
Sterling fe Witthoeft ( 2011 1 and Sterliug ( 2011). RR rate c oeffic ients for all other Z = 31- 
36 ions are extrapolations of the fits of Verne r & Ferland (11996 1 for Z < 30 species. The 
low-T DR rate coefficients for other Z > 30 species are treated as d escribed in iPaper I . 
High- T DR rate coef f icient s are taken from iMazzitelli fe Mattiolil (120021) for Ga and Ge and 
from Fournier et al. (120001 ) for Kr. The As, Se, and Br liigh-T DR rate coefficients are 
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assumed to be equal to the rate coefficient of the Kr ion with the same charge. While that 
approximation is crude, uncertainties in the high-T DR rate coefficients are unlikely to be 
important in models of photoionized nebulae, in which low-T DR is dominant. 


We utilize the Landau-Zener and Demkov CT rate coefficient calculations of Sterl ing fe Stanc il 


12 01111 for collisions between low-charge (up to five times ionized) Ge, Se, Br, and Kr ions 
with neutral hydrogen. CT between H + and neutral Ge, Se, Br, and Kr are also included. 
For the first four ions of Ga and As, the CT rate coefficients are taken to be the average of 
those for Cu and Ni ions of the same charge. For all other ionization states of Z = 31-36 
species, CT is treated analytically as described in Ferland et al. fjl9971 ) and Paper I. We do 
not include CT reactions with He, though in some cases these could affect the ionization 
balance. 


We assume that the Ga-Kr electron-impact ionization rate coefficients are equal to 
the corresponding ion of Zn (or H-like Zn for higher ionization states). Like high-T DR, 
collisional ionization does not play a significant role in photoionized plasmas such as PNe 
and H 11 regions. In addition, Auger ionization is disabled for all six elements. 

Ga-Kr have very low cosmic abundances, and should not affect the ionization balance 
or spectra of lighter elements. We computed models of the individual PNe in )J3] with the 
unaltered version of Cloudy to verify that the light element spectra are unaffected. We 
tested the implementation of Ga-Kr in the code by running the “smoke test” routine, which 
showed that the ionization distributions of these species are consistent with those for Z < 30 
elements. Finally, we ran the full test suite of routines to ensure that our modifications of 
Cloudy did not affect the code’s functionality for lighter species. 


3. MODELS OF INDIVIDUAL NEBULAE 


Our first objective is to test whether simulations using the new atomic data reproduce 
the observed ionization balance of Kr in PNe. Specifically, we model the intensities of Kr 
emission lines and compare them to the observed intensities of PNe in which multiple Kr 
ions have been detected. 


We have selected 15 PNe that exhibit emission lines from multiple Kr ions in their 
optical and near-infrared spectra (Tabled]). Ten of these obje cts were modeled in [Paper 1 . 
We draw the other objec ts from optical spectra published s ince IPaper II (ISharpee et al.ll2007 : 
Garcfa-Roias et al1l2012 1. with the exception of NGC 5315 ( Peimbert et ali2004 h Intensities 


for the K band emission lines from Paper II were determined using the extinction coefficients 
listed in Table Q] and the extinction law adopted by the optical data reference. We corrected 
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for aperture effects (e.g., different slit wi dths) by forcing the H I Bry/H f3 ratio to equal 
the theoretical value (Storey & Hummer 1995 ) for the nebula’s temperature and density. 
This line ratio is relatively insensitive to temperature and density in the regime of physical 
conditions found in PNe and H If regions. 


We modeled these 15 PNe with our modified version of Cloudy C13 ( Ferland et al.1 20131). 
adjusting the input parameters to best reproduce the observed optical and (when available) 
UV and IR spectra. The methodology is described in detail in Paper I], and here we provide 
only a brief synopsis. 

Our primary goal in this exercise is to simulate the ionization balance of each nebula, and 
not necessarily to produce a model that is physically accurate in terms of other nebular and 
stellar properties. Given this focus, we have made a number of simplifying assumptions in our 
models. We assume that the nebulae are spherical and have uniform densities, with covering 
and filling factors of unity. Furthermore, we do not incorporate dust physics, aside from 
the depletion of refractory elements. This assumption likely aff ects the thermal balance and 


derived central star temperatures of the nebular models (e.g., Dopita & Sut herla n d 2000)). 


However, we show in §4 that including dust does not significantly affect the ionization balance 
of Se or Kr. Despite these simplifications, the observed spectra are satisfactorily reproduced. 

Several input parameters must be specified in Cloudy models, including the central star 
temperature, gravity and luminosity; the inner and outer radii of the nebula; the nebular 
hydrogen density; and abundances of elements. For the ten PNe modeled in [Paper I . we used 
those input parameters as the starting values for our new models. The initial in put par ame¬ 


ters fo r the other ob jects were taken fromlSharpee et al.l (120071) f or IC 2501. IPottasch et al. 


(120051) for IC 4 191. iBernard Salas et al.l (120021) for NGC 2440, IPottasch et al.l (120021) for 
NGC 5315, and Pottasch fe Bernard-Salas ( 2008 ) for NGC 6369 . We obtained in i tial st ellar 
temperature and gravity data for NGC 5315 and NGC 2440 from Marcol ino et al. (120071 ) and 
Heapl (|1987l), respectively. Aside from stellar gravities and nebular inner radii, which were 
adjusted manually, all other modeled parameters were determined via the subplex optimiza¬ 
tion routine in Cloudy. In some cases, parameters were varied manually after optimization 
to improve the fit to the observed spectra. 

For the spectral energy distribution of t he PN central stars, we utilized NLTE stellar 
atmospheres including metal line blanketing (IRauchl 2 003) . Since thes e 15 P Ne are Galactic 
disk objects, we assumed solar metallicities for the central stars. The Rauch (120 030 grid of 
models spans the temperature range 50,000-190,000 K, and hence could not be used to model 
the low-ionization PN IC 418. For that object, we interpolated on two stellar atmosphere 
models ( Paper J ) computed by T. Rauch (2005, private communication) with 7j.fr = 30,000 K 
and 40,000 K and log g = 3.4. Uncertainties in T e s were estimated by manually adjusting 

























































Tabic 1. 


Observational Data References 


Object 

Name 

Optical 

Data Ref. 

c(m 

Other 

Ref. a 

IC 418 

Sharnee et al. 120031 

0.34 

Pottasch et al. (2004) 

IC 2501 

Sharnee et al. 120071 

0.55 

Goharii & Adams 119841: Milineo et al. 120021 

IC 4191 

Sharpee et al. 120071 

0.77 

Pottasch et al. 120051 

IC 5117 

Hyung et al. 120011 

1.40 


NGC 2440 

Sharnee et al. 120071 

0.55 

Bernard Salas et al. (2002) 

NGC 5315 

Pcimbert et al. (2004) 

0.74 

Pottasch et al. (2002) 

NGC 6369 

Garcia-Roias et al. (20121 

1.93 

Pottasch & Bernard-Salas 120081 

NGC 6572 

Liu et al. (2004a,bl 

0.48 

Hvune et al. 11994al b 

NGC 6741 

Liu et al. I2004a,bl 

1.15 


NGC 6790 

Liu et al. (2004a,bl 

1.10 

Aller et al. !l996l b 

NGC 6826 

Liu et al. (2004a.bl 

0.06 


NGC 6884 

Liu et al. I2004a,bl 

1.00 


NGC 6886 

Hvune et al. 119951 

0.90 

Pottasch & Surendiranath (20051 

NGC 7027 

Zhang et al. (2005) 

1.37 


NGC 7662 

Liu et al. (2004b,a) 

0.18 



Note. — Observational data references and extinction coeffici e nts c (H/3) are given. Near-infrared 
data are taken from Paper II . except for IC 5117 ( Rudv et al. 200 lh . NGC 7027, and NGC 7662 
( Geballe et al. 1991 b 

a Reference for abundances, UV, and mid-IR data if not from optical data reference. 

b [Kr IV] A5867.7 is taken from this reference and scaled to [Kr IV] A5346.0 of Liu et al. ( 2004al lbh. 
All other abundance, optical, and IR data (except the K band) are from Liu et al. ( 2004al fbl. 
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the central star temperature in each model, until diagnostic ratios and/or intensities of 
high-ionization species were no longer reproduced to within ~50-100%. 

The optimized input parameters to our models are given in Tables [2] and EJ Comparisons 
of modeled and observed line diagnostics and intensities are given in Tables HI [5j [6J and [7] 
with the exception of Kr lines which will be discussed in the following subsection. Line 
intensities are given on a normalized scale with J(H/3) = 100. For the majority of lines, 
the predicted intensities from the Cloudy models are within 25% of the observed intensities. 
Larger discrepancies exist for a handful of lines in some models, which are likely due to the 
simplifying assumptions noted above. 


Tabic 2. Model Input Parameters 


Parameter 

IC 418 

IC 2501 

IC 4191 

IC 5117 

NGC 2440 

NGC 5315 

NGC 6369 

NGC 6572 

Solar a 

T eff (10 3 K) 

38.8±3.0 

50.3±8.0 

118.8±10.0 

125.3±10.0 

188.4±30.0 

51.0±7.0 

62.0±10.0 

80.8±7.0 


log g 

3.4 

5.5 

6.0 

8.0 

6.7 

5.0 

5.1 

5.5 


L/Lq 

2330 

3850 

1680 

2610 

920 

990 

2620 

1480 


n n (10 3 cm -3 ) 

8.5 

11.6 

9.8 

59.0 

6.05 

30.2 

3.2 

19.2 


log(-R in /pc) 

-2.33 

-2.7 

-3.0 

-2.1 

-1.823 

-2.0 

-2.0 

-2.0 


log(_R out /pc) 

-1.13 

-1.405 

-1.524 

-1.68 

-1.324 

-1.234 

-1.18 

-1.44 


He/H 

0.105 

0.104 

0.112 

0.094 

0.124 

0.129 

0.101 

0.104 

0.085 

C/H 

8.67 

8.88 

8.78 

8.62 

8.52 

8.39 

8.20 

8.75 

8.43 

N/H 

7.92 

8.02 

8.22 

7.83 

8.65 

8.43 

8.03 

7.91 

7.83 

O/H 

8.58 

8.68 

9.09 

8.39 

8.64 

8.78 

8.65 

8.43 

8.69 

Ne/H 

7.90 

7.97 

8.40 

7.70 

7.93 

8.14 

7.95 

7.70 

7.93 

S/H 

6.63 

6.85 

7.24 

6.77 

6.59 

7.16 

6.84 

6.44 

7.12 

Cl/H 

4.29 

5.08 

5.35 

5.11 

5.05 

5.33 

4.94 

4.87 

5.50 

Ar/H 

6.13 

6.60 

6.46 

6.10 

6.45 

6.43 

6.24 

6.06 

6.40 

Se/H 




3.52±0.05 



3.66±0.14 

3.01±0.03 

3.34 

Kr/H 

3.93±0.23 

3.12±0.13 

3.22±0.30 

3.74±0.10 

3.44±0.10 

3.42±0.28 

3.94±0.18 

3.51±0.15 

3.25 


Note. — Input parameters for Cloudy models. All models were run with Cloudy version C13, modified as described in the text. The 
stellar effective temperature, gravity, and luminosity (relative to solar), nebular hydrogen density, and inner and outer nebular radii are 
given in the first six rows. Elemental abundances follow, with He/H on a linear scale, and the other elements given as 12+ log(X/H). 

a From IAsp lu ncTet all (120091) 
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Table 3. Model Input Parameters 


Parameter 

NGC 6741 

NGC 6790 

NGC 6826 

NGC 6884 

NGC 6886 

NGC 7027 

NGC 7662 

Solar a 

Teff (10 3 K) 

144.4+20.0 

96.7+10.0 

50.1+5.0 

114.0+10.0 

164.6+20.0 

162.7+20.0 

126.4+10.0 


log 9 

5.84 

5.5 

5.0 

5.7 

7.0 

6.0 

5.7 


L/Lq 

520 

5800 

12,400 

990 

780 

2930 

1550 


n h (10 3 cm -3 ) 

5.5 

38.5 

1.6 

7.0 

6.1 

55.0 

3.2 


log(i?in/pc) 

-2.24 

-2.3 

-1.60 

-2.7 

-2.79 

-2.05 

-1.60 


log(l? out /pc) 

-1.16 

-1.83 

-0.82 

-1.52 

-1.47 

-1.68 

-1.523 


He/H 

0.105 

0.101 

0.099 

0.099 

0.107 

0.100 

0.097 

0.085 

C/H 

8.57 

8.37 

8.29 

8.80 

8.87 

9.01 

8.59 

8.43 

N/H 

8.20 

7.57 

7.85 

8.23 

8.25 

8.23 

7.88 

7.83 

O/H 

8.58 

8.39 

8.64 

8.69 

8.66 

8.81 

8.57 

8.69 

Ne/H 

7.90 

7.73 

7.86 

7.92 

8.05 

8.04 

7.80 

7.93 

S/H 

6.77 

6.35 

6.42 

6.87 

6.88 

7.10 

6.76 

7.12 

Cl/H 

4.96 

4.84 

4.96 

5.12 

5.13 

5.33 

5.19 

5.50 

Ar/H 

6.32 

5.78 

6.11 

6.26 

6.37 

6.36 

6.11 

6.40 

Se/H 

3.31+0.11 

2.91+0.05 

2.95+0.13 

3.40+0.06 

3.52+0.08 

3.68+0.08 

3.53+0.10 

3.34 

Kr/H 

3.69+0.20 

3.26+0.10 

3.63+0.11 

3.69+0.15 

3.70+0.20 

4.14+0.11 

3.74+0.10 

3.25 


Note. — Input parameters for Cloudy models. All models were run with Cloudy version C13, modified as described 
in the text. The stellar effective temperature, gravity, and luminosity (relative to solar), nebular hydrogen density, and 
inner and outer nebular radii are given in the first six rows. Elemental abundances follow, with He/H on a linear scale, 
and the other elements given as 12+ log(X/H). 


^From Asplund et al. ( 2009h 
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Tablc 4. Predicted vs. Observed Diagnostics and Line Intensities 


Ion A 

IC 418 

IC 2501 

IC 4191 

IC 5117 

Pred. 

Obs. 

Pred. 

Obs. 

Pred. 

Obs. 

Pred. 

Obs. 

Ar Illna 

2.42E+02 

2.71E+02 

1.50E+02 a 

1.93E+02 a 

1.49E+02 a 

8.54E+01 a 

7.98E+01 

1.35E+02 

Ar IV nn 

1.61E+00 

1.20E+00 

1.87E+00 

1.50E+00 

1.73E+00 

1.68E+00 

4.16E+00 

4.18E+00 

Cl Illnn 

6.07E-01 

5.10E-01 

5.29E-01 

5.61E-01 

5.67E-01 

4.92E-01 

2.95E-01 

3.28E-01 

N Ilna 

8.41E+01 

7.84E+01 

5.68E+01 

5.64E+01 

6.26E+01 

4.29E+01 

1.81E+01 

2.53E+01 

O Ilna 

8.54E+00 

6.98E+00 

5.30E+00 

4.36E+00 

6.22E+00 

2.65E+00 

9.49E-01 

6.95E-01 

O linn 

2.18E+00 

2.37E+00 

2.26E+00 

2.38E+00 

2.23E+00 

2.31E+00 

2.55E+00 

3.10E+00 

O Ulna 

3.19E+02 

3.08E+02 

2.45E+02 

2.39E+02 

2.34E+02 

2.13E+02 

7.16E+01 

8.83E+01 

S linn 

4.98E-01 

4.70E-01 

4.96E-01 

5.00E-01 

4.81E-01 

5.13E-01 

4.46E-01 

4.36E-01 


Note. — Table [4] is published in its entirety in the electronic edition of the article. A portion is shown here for guidance 
regarding its form and content. Comparison of predicted and observed line intensities (on the scale /(H/3) = 100) and 
diagnostic ratios. Wavelengths are in A except where noted. The diagnostic ratios are defined as follows: Ar Illna = [Ar mi 
(7135 + 7751)/5192; Ar IVnn = [Ar IV] 4740/4711; Cl Illnn = [Cl III] 5518/5538; N Una = [N II] (6548 + 6584)/5755; 
O Ilna = [O II] (3726 + 3729)/(7320 + 7330); O linn = [O II] 3726/3729; O Ulna = [O III] (4959 + 5007)/4363; and 
S linn = [S II] 6716/6731. The notation “na” indicates a nebular to auroral line ratio, and “nn” indicates nebular to 
nebular. Nebular lines arise from the first excited term of the ground configuration whereas auroral lines come from the 
second excited term. 

a Ar I Ilna = [Ar III] 7135/5192. 
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Table 5. Predicted vs. Observed Diagnostics and Line Intensities 


Ion 

A 

NGC 2440 

NGC 5315 

NGC 6369 

NGC 6572 

Pred. 

Obs. 

Pred. 

Obs. 

Pred. 

Obs. 

Pred. 

Obs. 

ArIIlna 


7.60E+01 a 

8.48E+01 a 

2.15E+02 

3.00E+02 

1.52E+02 

2.18E+02 

9.93E+01 

1.69E+02 

Ar IV nn 


1.34E+00 

1.16E+00 

3.22E+00 

4.69E+00 

1.07E+00 

1.19E+00 

2.35E+00 

1.98E+00 

Cl TTTnn 


7.16E-01 

7.50E-01 

3.52E-01 

3.41E-01 

8.94E-01 

8.24E-01 

4.37E-01 

4.40E-01 

N Una 


5.66E+01 

4.67E+01 

4.24E+01 

4.50E+01 

6.58E+01 

4.34E+01 

3.58E+01 

3.78E+01 

O Una 


7.98E+00 

8.39E+00 

2.40E+00 

2.21E+00 

1.32E+01 

5.55E+00 

2.80E+00 

3.06E+00 

O Unn 


1.99E+00 

1.79E+00 

2.49E+00 

2.67E+00 

1.74E+00 

1.81E+00 

2.37E+00 

2.48E+00 

O Ulna 


8.35E+01 

7.14E+01 

2.72E+02 

2.56E+02 

1.90E+02 

1.64E+02 

2.70E+01 b 

4.30E+01 b 

S Unn 


5.60E-01 

6.45E-01 

4.58E-01 

4.76E-01 

5.75E-01 

5.98E-01 

4.76E-01 

4.70E-01 


Note. — Table [5] is published in its entirety in the electronic edition of the article. A portion is shown here for guidance 
regarding its form and content. Comparison of predicted and observed line intensities (on the scale /(H/3) = 100) and 
diagnostic ratios. Wavelengths are in A except where noted. The diagnostic ratios are defined as follows: Ar Illna = [Ar iii] 
(7135 + 7751)/5192; Ar IVnn = [Ar IV] 4740/4711; Cl TT T nn = [Cl III] 5518/5538; N Una = [N II] (6548 + 6584)/5755; 
O Ilna = [O II] (3726 + 3729)/(7320 + 7330); O linn = [O II] 3726/3729; O Ulna = [O III] (4959 + 5007)/4363; and 
S linn = [S II] 6716/6731. The notation “na” indicates a nebular to auroral line ratio, and “nn” indicates nebular to 
nebular. Nebular lines arise from the first excited term of the ground configuration whereas auroral lines come from the 
second excited term. 

a Ar Ulna = [Ar III] 7135/5192 

b O IHna = [O III] 4959/4363 
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Tablc 6. Predicted vs. Observed Diagnostics and Line Intensities 


Ion 

A 

NGC 6741 

NGC 6790 

NGC 6826 

NGC 6884 

Pred. 

Obs. 

Pred. 

Obs. 

Pred. 

Obs. 

Pred. 

Obs. 

ArIIlna 


9.71E+01 

1.32E+02 

7.68E+01 

9.16E+01 

2.02E+02 

2.56E+02 

1.19E+02 

1.51E+02 

Ar IV nn 


1.28E+00 

1.26E+00 

3.35E+00 

4.06E+00 

9.04E-01 

9.32E-01 

1.43E+00 

1.56E+00 

ClIllnn 


7.53E-01 

7.30E-01 

3.37E-01 

4.07E-01 

1.08E+00 

1.10E+00 

6.69E-01 

7.04E-01 

N Ilna 


5.44E+01 

6.48E+01 

2.29E+01 

2.27E+01 

8.68E+01 

7.91E+01 

5.03E+01 

5.25E+01 

O ilna 


8.46E+00 

9.47E+00 

1.38E+00 

1.51E+00 

2.56E+01 

1.48E+01 

6.74E+00 

7.02E+00 

O linn 


1.93E+00 

2.04E+00 

2.50E+00 

2.45E+00 

1.43E+00 

1.48E+00 

2.09E+00 

2.17E+00 

O Ulna 


2.33E+01 a 

2.67E+01 a 

1.87E+01 a 

2.24E+01 a 

6.43E+01 a 

6.71E+01 a 

3.41E+01 a 

3.78E+01 a 

S linn 


5.62E-01 

5.70E-01 

4.55E-01 

4.57E-01 

6.74E-01 

7.30E-01 

5.02E-01 

5.32E-01 


Note. — Table [6] is published in its entirety in the electronic edition of the article. A portion is shown here for guidance 
regarding its form and content. Comparison of predicted and observed line intensities (on the scale /(H/3) = 100) and 
diagnostic ratios. Wavelengths are in A except where noted. The diagnostic ratios are defined as follows: Ar Illna = [Ar mi 
(7135 + 7751)/5192; Ar IVnn = [Ar IV] 4740/4711; Cl Illnn = [Cl III] 5518/5538; N Una = [N II] (6548 + 6584)/5755; 
O Ilna = [O II] (3726 + 3729)/(7320 + 7330); O linn = [O II] 3726/3729; O Ulna = [O III] (4959 + 5007)/4363; and 
S linn = [S II] 6716/6731. The notation “na” indicates a nebular to auroral line ratio, and “nn” indicates nebular to 
nebular. Nebular lines arise from the first excited term of the ground configuration whereas auroral lines come from the 
second excited term. 

a O Ulna = [O III] 5007/4363. 
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Tabic 7. Predicted vs. Observed Diagnostics and Line Intensities 


Ion 

A 

NGC 6886 

NGC 7027 

NGC 7662 

Pred. 

Obs. 

Pred. 

Obs. 

Pred. 

Obs. 

Ar IHna 


1.04E+02 

1.09E+02 

1.06E+02 

9.79E+01 

9.97E+01 

9.68E+01 

Ar IV nn 


1.35E+00 

1.29E+00 

4.12E+00 

3.70E+00 

1.06E+00 

1.26E+00 

Cl Illnn 


7.11E-01 

3.94E-01 

2.99E-01 

2.90E-01 

8.97E-01 

9.87E-01 

N Ilna 


5.22E+01 

6.38E+01 

2.26E+01 

2.20E+01 

5.05E+01 

5.05E+01 

O Ilna 


7.50E+00 

4.50E+00 

1.16E+00 

8.70E-01 

1.14E+01 

8.63E+00 

O linn 


2.02E+00 

1.53E+00 

2.55E+00 

2.76E+00 

1.73E+00 

1.70E+00 

O Ulna 


1.04E+02 

1.22E+02 

9.11E+01 

8.75E+01 

2.59E+01 a 

2.33E+01 a 

S linn 


5.15E-01 

4.84E-01 

4.47E-01 

4.40E-01 

5.79E-01 

6.58E-01 


Note. — Table [7] is published in its entirety in the electronic edition of the article. A portion is 
shown here for guidance regarding its form and content. Comparison of predicted and observed line 
intensities (on the scale /(H/3) = 100) and diagnostic ratios. Wavelengths are in A except where 
noted. The diagnostic ratios are defined as follows: Ar IHna = [Ar iii] (7135 + 7751)/5192; 
Ar IVnn = [Ar IV] 4740/4711; Cl Illnn = [Cl III] 5518/5538; N Una = [N II] (6548 + 
6584)/5755; O Una = [O II] (3726 + 3729)/(7320 + 7330); O linn = [O II] 3726/3729; O IHna 
= [O III] (4959 + 5007)/4363; and S linn = [S II] 6716/6731. The notation “na” indicates a 
nebular to auroral line ratio, and “nn” indicates nebular to nebular. Nebular lines arise from the 
first excited term of the ground configuration whereas auroral lines come from the second excited 
term. 

a O IHna = [O III] 4959/4363. 
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We determined uncertainties in the modeled Se and Kr abundances (Tables [2] and [3]) by 
considering errors in line intensities and central star temperature. For the upper and lower 
limits to the stellar temperature, we re-optimized the Se and Kr abundances to estimate their 
uncertainty. The effects of uncertainties in line intensities on the abundance determinations 
were treated differently for Se and Kr. In the case of Kr, we manually adjusted each model’s 
Kr abundance until the most discrepant Kr line intensity was reproduced. We adopted the 
difference with the best model as the uncertainty, or 25% (the assumed error bars in the 
optical line intensities), whichever value was larger. In contrast, only one line of Se was 
detected (if at all) in the modeled nebulae, and we utilized the percentage uncertainty in the 
intensity as that of the abundance (this underestimates the magnitude of the actual error 
bars, as discussed in Paper I). The overall abundance uncertainties were calculated by adding 
those from the intensities and variations in the central star temperature in quadrature. 


3.1. The Ionization Balance of Krypton 


In Table [8J, we compare the observed Kr line intensities with those predicted by our 
models. We find that the ionization balance of Kr is not reproduced well by our models that 
incorporate the new atomic data. Specifically, the models predict values for the ionic ratio 
Kr 3+ /Kr 2+ that are lower than the observed values. [Kr III] line strengths are predicted to be 
larger than observed by about 50% up to factors of 2-3, while the [Kr IV] lines are predicted 
to be weaker than observed by similar amounts. This is a systematic effect, in that it is 
seen in all of the 15 modeled PNe independent of ionization level, central star temperature, 
morphology, extinction, density, and other nebular properties. This result suggests that the 
errors in the modeled Kr ionic ratios are due to errors in the new atomic data. To verify 
this, we first investigate whether observational uncertainties or modeling biases can explain 
the poorly-reproduced Kr ionization equilibrium. 


The fluxes of weak, low signal-to-noise emission lines are often overestimated compared 
to stronger lines (e.g., Rola & Pelat 1994). Because of the low cosmic abundance of Kr, 
even “strong” Kr emission lines are comparable in flux to weak permitted features of more 
abundant elements. As seen in Table 0 [Kr III] and [Kr V] lines tend to be weaker than 
[Kr IV] lines for PNe of moderate to high excitation. Therefore, if line measurement errors 
substantially affect the Kr ionization equilibrium results, one would expect the observed 
[Kr III] and [Kr V] lines to be stronger than predicted by models. However, the opposite 
effect is seen in our models, where the predicted [Kr III] line intensities are too large compared 
to observational measurements. In NGC 7027, the observed [Kr V] intensities are also much 
weaker than predicted. Therefore bias in the measurement of weak emission lines does not 






Table 8. Predicted vs. Observed Krypton Line Intensities 


Ion 

A 

Pred. a 

Pred. b 

Obs. 

Pred. a 

Pred. b 

Obs. 

Pred. a 

Pred. b 

Obs. 




IC 418 



IC 2501 



IC 4191 


[Kr ill] 

6826 

4.71E-02 

2.83E-02 

2.90E-02 

8.08E-03 

5.15E-03 

4.60E-03 

6.38E-03 

4.60E-03 

1.50E-03 

[Kr III] 

2.199/xm 

1.44E-01 

8.75E-02 

6.99E-02 

2.19E-02 

1.40E-02 


1.72E-02 

1.23E-02 


[Kr IV] 

5346 

8.55E-04 

1.88E-03 

3.50E-03 

4.04E-03 

5.54E-03 

6.10E-03 

9.78E-03 

1.24E-02 

2.40E-02 

[Kr IV] 

5868 

1.30E-03 

2.85E-03 

3.50E-03 

6.08E-03 

8.34E-03 

8.50E-03 

1.43E-02 

1.81E-02 

1.90E-02 

[Kr V] 

6256 

1.54E-10 

1.74E-12 


1.31E-08 

4.35E-10 


5.40E-04 

1.03E-04 


[Kr V] 

8243 

1.40E-10 

1.59E-12 


1.19E-08 

3.97E-10 


4.93E-04 

9.37E-05 


[Kr VI] 

1.233/xm 




8.03E-17 



9.96E-05 

3.72E-05 





IC 5117 



NGC 2440 



NGC 5315 


[Kr III] 

6826 

4.13E-02 

3.05E-02 


4.20E-02 

1.03E-02 


1.75E-02 

1.17E-02 

6.20E-03 

[Kr III] 

2.199/im 

7.16E-02 

5.26E-02 

4.23E-02 

8.82E-02 

2.08E-02 


5.08E-02 

3.44E-02 


[Kr IV] 

5346 

8.76E-02 

1.07E-01 

1.29E-01 

4.65E-02 

3.48E-02 

3.50E-02 

3.33E-03 

4.99E-03 

9.40E-03 

[Kr IV] 

5868 

1.23E-01 

1.49E-01 

1.45E-01 

6.25E-02 

4.62E-02 

4.60E-02 

5.19E-03 

7.78E-03 


[Kr V] 

6256 

1.75E-03 

2.20E-04 


2.78E-02 

4.77E-03 

4.40E-02 

1.48E-09 

9.74E-14 


[Kr V] 

8243 

1.59E-03 

2.01E-04 


2.53E-02 

4.35E-03 


1.35E-09 

8.88E-14 


[Kr VI] 

1.233/xm 

5.24E-05 

6.84E-06 


5.22E-02 

1.18E-02 








NGC 6369 



NGC 6572 



NGC 6741 


[Kr III] 

6826 

5.58E-02 

3.92E-02 

2.10E-02 

2.60E-02 

1.87E-02 

1.80E-02 

4.20E-02 

2.45E-02 

3.70E-02 

[Kr III] 

2.199/rm 

1.36E-01 

9.50E-02 


5.07E-02 

3.62E-02 

2.49E-02 

8.37E-02 

4.81E-02 

5.40E-02 

[Kr IV] 

5346 

4.20E-02 

5.68E-02 

8.80E-02 

2.97E-02 

3.90E-02 

4.70E-02 

4.23E-02 

6.61E-02 

6.90E-02 

[Kr IV] 

5868 

6.08E-02 

8.22E-02 

1.23E-01 

4.14E-02 

5.44E-02 

6.40E-02 

5.68E-02 

8.81E-02 

5.50E-02 

[Kr V] 

6256 

2.52E-06 

1.49E-07 


8.69E-06 

5.85E-07 


1.81E-02 

5.71E-03 


[Kr V] 

8243 

2.30E-06 

1.36E-07 


7.92E-06 

5.34E-07 


1.65E-02 

5.20E-03 


[Kr VI] 

1.233/im 

5.14E-04 

2.82E-15 


6.34E-13 

4.08E-14 


1.08E-02 

5.10E-03 





NGC 6790 



NGC 6826 



NGC 6884 


[Kr III] 

6826 

1.34E-02 

9.91E-03 

1.20E-02 

2.67E-02 

1.57E-02 

1.60E-02 

2.76E-02 

1.99E-02 

1.50E-02 

[Kr III] 

2.199jum 

1.96E-02 

1.69E-02 


7.51E-02 

4.41E-02 


5.95E-02 

4.26E-02 


[Kr IV] 

5346 

3.33E-02 

4.70E-02 

5.30E-02 

1.80E-02 

2.38E-02 

2.30E-02 

5.37E-02 

6.58E-02 

9.50E-02 

[Kr IV] 

5868 

4.58E-02 

6.47E-02 

5.20E-02 

2.66E-02 

3.50E-02 

4.50E-02 

7.45E-02 

9.13E-02 

1.14E-01 

[Kr V] 

6256 

1.94E-04 

1.79E-05 


5.15E-08 

1.20E-09 


1.39E-03 

1.97E-04 


[Kr V] 

8243 

1.77E-04 

1.63E-05 


4.70E-08 

1.09E-09 


1.26E-03 

1.80E-04 


[Kr VI] 

1.233/xm 

2.75E-06 

2.85E-07 


4.36E-17 



1.40E-04 

3.48E-05 





NGC 6886 



NGC 7027 



NGC 7662 


[Kr III] 

6826 

4.34E-02 

2.10E-02 


7.59E-02 

4.59E-02 

4.10E-02 

2.06E-02 

1.25E-02 

1.30E-02 

[Kr III] 

2.199/rm 

8.89E-02 

4.28E-02 

5.89E-02 

1.53E-01 

9.20E-02 

8.74E-02 

4.18E-02 

2.54E-02 


[Kr IV] 

5346 

4.88E-02 

6.25E-02 

8.00E-02 

9.05E-02 

1.48E-01 

1.48E-01 

1.07E-01 

1.30E-01 

1.59E-01 

[Kr IV] 

5868 

6.64E-02 

8.43E-02 

5.40E-02 

1.30E-01 

2.10E-01 

2.22E-01 

1.45E-01 

1.75E-01 

1.60E-01 

[Kr V] 

6256 

2.45E-02 

6.87E-03 


4.47E-02 

2. HE—02 

2.50E-02 

1.61E-02 

2.44E-03 


[Kr V] 

8243 

2.24E-02 

6.27E-03 


4.08E-02 

1.93E-02 

1.40E-02 

1.47E-02 

2.23E-03 


[Kr VI] 

1.233/xm 

2.64E-02 

1.10E-02 


1.11E-01 

6.50E-02 


1.42E-03 

2.37E-04 



Note. — Comparison of predicted and observed Kr line intensities (on the scale 7(H/3) = 100). Wavelengths are in A unless noted. 
a Predicted Kr line intensities, using the photoionization and recombination data from ISterlind d201lh . 

b Predicted Kr line intensities, after adjusting the PI cross sections and DR rate coefficients within their uncertainties (see text). 
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explain the inconsistency between the predicted and observed Kr ionization balance. 


Calibration uncertainties can also potentially be important, particularly since many of 
the optical data we utilize are cross-dispersed echclle spectra, for which flux calibrations are 
often difficult and uncertain. However, since the data we utilize were taken with different 
instruments and reduced by independent teams of researchers, one would expect that calibra¬ 
tion errors would be random and should not produce systematic effects. Even if systematic 
effects were introduced by improper calibrations, then these should in principle be seen for 
other species. However, no such trends are seen for other elements. 


Assumptions and/or uncertainties in our numerical models are another potential expla¬ 
nation. Indeed, we make several simplifying assumptions in our models: we assume that 
the nebulae are uniform density and isotropic, we ignore the effects of dust on the ionized 
gas, and we assume that the filling and covering factors are equal to unity. These assump¬ 
tions certainly affect the predicted emission line spectra to some extent. Discrepancies in 
the predicted S, Ar, and Cl ionization balance for some objects are likely due in part to 
our simplifications. However, it is important to note that there are no systematic trends 
seen in the predicted ionization balance of elements other than Kr. This suggests that the 
simplifying assumptions in our models are unlikely to be the source of the Kr ionization 
equilibrium disparity. The fact that the inclusion or exclusion of dust in our models does 
not substantially affect the Kr ionization balance (j|4]) further supports this argument. 


Therefore it appears to be unlikely that observational or calibration errors and simpli¬ 
fications in our models can explain the disagreement between the predicted and observed 
Kr ionization balance. On the other hand, it is possible that errors and uncertainties in 
the atomic data may be responsible for the ionization balance discrepancy. ISterlina (2011) 
cite error bars in their computed PI cross sections and RR and DR rate coefficients, based 
on comparisons to experimental measurements, the electronic configurations included in the 
calculations, and the effects of variations in the radial scaling parameters of the orbital wave- 
functions. For example, the computed Kr 3+ direct PI cross sec tion is a factor of two l arge r in 
magn itude than the experimental PI cross section of Lu et al. (2006a) (see Fig. 5 of [Sterling 


2 0111) . and efforts to resolve this discrepancy were not successful. For other ions, the PI cross 
sections were cited to be accurate to within 30-50%. DR rate coefficients for near-neutral 
heavy element ions have large uncertainties regardless of the theoretical method utilized, 
due to the unknown energies of autoionizing resonances near the ionization threshold. In 
the “low” temperature environments of photoionized nebulae, electrons are captured into 
these levels, and hence whether they lie just above the ionization threshold or just below 
can substantia lly affect the computed DR rate coefficient flSavin et al.l Il999l: iFerlandl 120031) . 


Sterling 12 011 ) estimated uncertainties of factors of 2-3 for DR of Kr ions, though the actual 
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errors could be larger. 

Uncertainties in DR rate coefficients are likely to be more important than those for 
RR rate coefficients, since DR is the dominant recombination mechanism for Kr ions at 
electron temperatures near 10 4 K ([Sterling! 2 01 It ) . Moreover, the RR rate coefficients have a 
relatively small uncertainty (< 10 % for all but singly-ionized Kr) compared to DR. The CT 
calculations of [Sterling fe Stancill f 201 ll) are another potential source of error, as Landau- 
Zener calculations are only accurate to within a factor of ~3 for systems with large CT rate 
coefficients. However, the importance of CT decreases with increasing charge of the reactant 
ion, as less neutral hydrogen is generally intermixed with more highly-charged species. 


Therefore, we consider the DR and PI data to be the most potentially important sources 
of uncertainty in the Kr ionization balance. In Cloudy, we manually adjusted the PI cross 
sections and DR rate coefficients for various Kr ions to investigate how they affect the relative 
strengths of [Kr III] and [Kr IV] emission lines. We found that altering only the PI cross 
sections (by up to a factor of two) or only the DR rate coefficients (by up to a factor of 10) did 
not eliminate the systematic discrepancy in the predicted Kr ionization balance. However, by 
adjusting both the PI and DR data for Kr + , Kr 2+ , and Kr 3+ by modest amounts, systematic 
trends in the relative strengths of [Kr III] and [Kr IV] line intensities could be eliminated. 


The PI cross sections were altered within the uncertainties cited by Sterling (2011), while 
the DR rate coefficients were modified by an order of magnitude (larger than the cited - 
but approximate - uncertainties; see Table®. Specifically, we adjusted the PI cross sections 
in the following manner: Kr + was reduced by 50%, Kr 2 " 1 " was increased by 50%, and Kr 3+ 
was reduced by a factor of two. The larger ch ange to the Kr 3+ PI cross section is motivated 
by the experimental cross section (L u et a l. 2 006b ). which has a magnitude half that of the 
cross section computed by [Sterling] (20111). The DR rate coefficient of Kr 3+ (forming Kr 2+ ) 
was reduced by a factor of ten, while those of Kr 2+ and Kr 4+ were increased by a factor of 10. 
These changes reduced the Kr 2+ ionic fractions in our numerical simulations, while increasing 
those of Kr 3+ . The alterations to the atomic data also reduced the (initially overestimated) 
[Kr V] AA6256, 8243 line strengths in NGC 7027, though the agreement worsened for A6256 
in NGC 24400 


With these changes to the Kr PI and DR atomic data, there is no systematic disagree¬ 
ment between the modeled and observed Kr ionization balance. The remaining inconsisten¬ 
cies in the predicted Kr line intensities are similar in magnitude to those seen for S and Ar, 
and can be attributed to observational uncertainties and/or model simplifications. 


1 Sharpee et al. (j2007i) note that this line is affected by telluric contamination, and hence its flux may be 


uncertain. 
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Table 9. 


Adjustments made to Kr atomic data of Sterling (2011) 


Ion 

<Tpi (adjusted)/ 

CTpi(Sll) 

«dr (adjusted)/ 
«dr(S11) 

Kr+ 

0.67 


Kr 2+ 

1.5 

10.0 

Kr 3+ 

0.5 

0.1 

Kr 4+ 


10.0 


Note. — The adjustm ent factors mad e 
to the Kr atomic data of I Ster ling (2011). 
The first column lists the ion before the 
indicated reaction occurs. Therefore, the 
PI cross section (cr pi ) of Kr + refers to pho- 
toionization of Kr + to form Kr 2+ . Like¬ 
wise, the DR rate coefficient (ckdr) for 
Kr 2+ refers to Kr 2+ forming Kr + . 
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This exercise highlights the value of testing theoretical calculations against experimental 
and observational measurements, particularly for complex systems such as near-neutral ions 
of heavy elements. As discussed in the introduction, PNe have served as natural laboratories 
to evaluate the accuracy of atomic data determinations for many ions. Observations, paired 
with increasingly sophisticated theoretical and experimental methods to determine atomic 
properties, have led to more reliable atomic data for many species. Therefore it should come 
as little surprise that a similar scenario occurs for Kr. 


Because the atomic data governing the ionization balance of Se cannot currently be 
tested in the absence of detections of ions other than Se 3+ , derived Se abundances will have 
inherently larger uncertainties than Kr. Ho wever, s ince Se 3+ is expected to be the dominant 
Se ion in moderate to high-excitation PNe (IDinersteinl 200111 . the corrections for unobserved 
ions should be relatively modest and thus the effects of atomic data uncertainties mitigated. 


3.2. Comparison of Modeled Se and Kr Abundances with Previous Results 


We compare the derived Se and Kr abundances from our new Cloudy models with those 
of Paper I in Table flOl For Kr, we show the abundances determined both before and after we 
adjusted the PI and DR data to better reproduce the Kr ionization balance. The abundances 
found using the new atomic data show remarkably good agreement with those found using 
the approximate atomic data of Pap er I. 


To explore the reason for this, in Table [TT] we compare t he PI cross sec ti ons near the 
ioni zation threshol d and DR rate coefficients at 10 4 K from [Sterling fe Witthoeftl (1201111 
and iSterlingl (20.11) to the approximate values of Pa p er ij (we do not consider RR, which at 
typical PN temperatures is dominated by DR for these species). For both Se and Kr, the PI 
cross sections are generally smaller (with the exception of singly-ionized species) while the 
DR rate coefficients are generally larger (with the exception of doubly-ionized species). 

In the case of Kr, the DR rate c oefficients of Kr 3+ (forming Kr 2+ ) and Kr 4+ are sig¬ 
nificantly larger than used in Pap er Ij. Combined with the smaller PI cross sections of these 
species, it is expected that the overall Kr ionization balance will shift to lower charge states. 
In Table [T21 we compare the modeled Kr ionic fractions (averaged over volume) with those 
of Paper I for each PN common to both studies. In all models, the Kr ionic fractions are 
indeed shifted to less ionized species. The Kr 2+ fractional abundances are larger in our new 
models with the exception of IC 418, in which Kr + is now predi cted to be more abundant 
than Kr 2+ . Likewise, the Kr 3+ fraction is reduced compared to Paper I in all but the most 
highly-ionized nebulae, for which the Kr 3+ abundance is larger at the expense of Kr 4+ and/or 
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Table 10. Comparison of Modeled Se and Kr Abundances 


PN Name 

12 + log(Se/H) 


12 ± log(Kr/H) 


This work 

Paper I 

Adopted 

No Adjustment 

Paper I 

IC 418 



3.93T0.23 

3.95 

3.90T0.21 

IC 2501 



3.12T0.13 

3.12 


IC 4191 



3.22T0.30 

3.20 


IC 5117 

3.52T0.05 

3.67T0.05 

3.74T0.10 

3.72 

3.76T0.11 

NGC 2440 



3.44T0.10 

3.90 


NGC 5315 



3.42T0.28 

3.43 


NGC 6369 

3.66T0.14 


3.94T0.18 

3.92 


NGC 6572 

3.01T0.03 

3.11T0.08 

3.51T0.15 

3.47 

3.54T0.18 

NGC 6741 

3.31T0.11 

3.39T0.12 

3.69T0.20 

3.75 

3.75T0.13 

NGC 6790 

2.91T0.05 

3.10T0.05 

3.34T0.10 

3.23 

3.24T0.14 

NGC 6826 

2.95T0.13 

2.89T0.13 

3.63T0.11 

3.67 

3.57T0.20 

NGC 6884 

3.40T0.06 

3.45T0.10 

3.69T0.15 

3.67 

3.65T0.16 

NGC 6886 

3.70T0.08 

3.59T0.11 

3.70T0.20 

3.87 

3.81T0.18 

NGC 7027 

3.68 ±0.08 

3.78T0.12 

4.14T0.11 

4.23 

4.26T0.12 

NGC 7662 

3.53T0.10 

3.50T0.10 

3.74T0.10 

3.79 

3.84T0.14 


Note. — Comparison o f Se and Kr abundances derived from models of 15 individ¬ 
ual PNe with those from Pap er I. For Kr, we show the adopted Kr abundances (after 
the atomic data was altered to improve the modeled Kr ionization balance in each 
PN), and those found before the atomic data were adjusted. The adjusted atomic 
data do not significantly impact Kr abundances when multiple Kr ions are detected, 
but will likely have a larger effect when a single ionization stage is detected. 
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Table 11. 


Comparison of Se and Kr Atomic Data with Paper I 


Ion 

dpi (thresh)/ 
dpi (Paper I) 

a DR (10 4 K)/ 
ann fPaper I) 

Se+ 

2.5 

27.1 

Se 2+ 

0.9 

0.3 

Se 3+ 

0.3 

1.1 

Se 4+ 

0.3 

5.5 

Kr+ 

1.6 

1.6 

Kr 2+ 

0.6 

0.2 

Kr 3+ 

0.7 

7.9 

Kr 4+ 

0.3 

24.4 


Note. — Ratio of PI cross 
sections (just above the ionization 


of 

Sterling; & Witthoeft 

Sterling! ( 

2011 

) relative 


proximate data of Paper I. The first 
column lists the ion before the indi¬ 
cated reaction occurs. Therefore, the 
PI cross section (cr pi ) of Se + refers 
to photoionization of Se + to form 
Se 2+ . Likewise, the DR rate coeffi¬ 
cient («dr) for Se 2+ refers to Se 2+ 
forming Se + . 
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Kr 5+ in our new models. 


The excellent agreement between the Kr abundances in our new models and previous 
results can be understood given that we selected PNe in which multiple Kr ions have been 
detected. In most objects, the detected species are Kr 2+ and Kr 3+ . These ions have an 
ionization potential range of 24.4-50.9 eV, and hence are expected to be the dominant Kr 
ions in all but the lowest and highest excitation PNe. Indeed, according to our Cloudy models 
the combined ionic fractions of Kr 2+ and Kr 3+ exceed 70-80% of the total Kr elemental 
abundance in all but IC 418, NGC 2440, NGC 67 41, a nd NGC 7662. These values are within 
5-10% of the combined ionic fractions from Paper I fTable fl2l) . Therefore, the corrections 
for unobserved Kr ions are small due to the detection of the dominant ionization stages. 

The agreement is not expected to be as strong when only one ionization stage is de¬ 
tected, as is the case for most of the PNe in the sample of Paper II. Indeed, the larger Kr 2+ 
ionic fractions in our modeled objects suggest that previous ICF schemes systematically 
overestimated total elemental Kr abundances in PNe. In fj5l we see that this is indeed the 


case. 


It is also worth noting that our adjustment of the Kr atomic data has very little ef¬ 
fect on the derived Kr elemental abundance in our models, due to the detection of both 
Kr 2+ and Kr 3+ (Table flOl) . The difference in the derived Kr abundance is greater than 
0.1 dex only in the case of NGC 2440 and NGC 6886. For NGC 2440, we were unable to 
fit [Kr V] A6256.1 simultaneous ly with [Kr IV] AA5346.0, 5867.7 after the atomic data ad¬ 
justments. I Shar pee et ah (120071) note that A6256.1 is blended with telluric features in their 
NGC 2440 spectrum, and hence its flux may be uncertain. Regardless, a larger Kr abun¬ 
dance was required to reproduce the Kr 3+ emission lines in NGC 2440 after the adjustment 
to the atomic data. For NGC 6886, our adjustment of the atomic data shifted a significant 
amount of Kr 2+ and Kr 4+ to Kr 3+ in our models. Along with the uncertainties in the [Kr IV] 
line intensities (see Table [8]) this required a larger Kr abundance to reproduce the observed 
Kr emission lines. 


Our Se abundances also agree well with those derived in Paper I. At face value, this 
seems more surprising, given that only a single ion of Se was detected. In many of the 15 
modeled PNe in which it was detected, Se 3+ comprises at least 50% of the total gaseous Se 
abundance, with the exceptions of NGC 6741, NGC 6886, NGC 7027, and NGC 7662. But 
even in those four PNe, the modeled Se 3+ fractional abundances a gree w ith those of Paper I 
to within 6%. For other objects, the ionic fraction differs from the Paper I models by up to 
25%. 


However, the reason for this seeming coincidence emerges when the new PI and DR 
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Table 12. Comparison of Kr Ionic Fractions With 


PN 

Name 

Kr + 

Kr 2+ 

Kr 3 + 

Kr 4 + 

Kr 5 + 

This work 

Paper I 

This work 

Paper I 

This work 

Paper I 

This work 

Paper I 

This work 

Paper I 

IC 418 

0.56 

0.40 

0.41 

0.55 

0.02 

0.03 





IC 5117 

0.10 

0.12 

0.26 

0.17 

0.61 

0.68 


0.02 



NGC 6572 

0.16 

0.13 

0.32 

0.21 

0.49 

0.64 





NGC 6741 

0.25 

0.29 

0.25 

0.20 

0.35 

0.30 

0.06 

0.09 

0.03 

0.09 

NGC 6790 

0.08 

0.09 

0.21 

0.15 

0.69 

0.74 


0.02 



NGC 6826 

0.08 

0.01 

0.40 

0.31 

0.52 

0.68 





NGC 6884 

0.06 

0.01 

0.28 

0.19 

0.66 

0.74 


0.04 


0.02 

NGC 6886 

0.14 


0.26 

0.11 

0.43 

0.47 

0.09 

0.16 

0.08 

0.17 

NGC 7027 

0.11 

0.12 

0.20 

0.14 

0.39 

0.33 

0.09 

0.10 

0.16 

0.19 

NGC 7662 

0.02 


0.13 

0.10 

0.82 

0.77 

0.03 

0.09 


0.04 


Note. — Comparison of modeled Kr ionic fractions (averaged over volume) between this work and Paper I. Ionic fractions less 


than 1% are omitted. 


to 

cn 
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data for Se ions are compared to the previous approximate values. Note in Table |TT] that 
the PI cross section _of Se 2-f ~ and DR rate coefficient of Se 3+ forming Se 2+ are quite similar 
to those adopted in Paper I. On the other hand, the PI cross section of Se 3+ is considerably 
smaller than 


Paper j ’s while the DR rate coefficient of Se 4+ forming Se 3+ is much larger. 


Therefore, it may be expected that the Se abundances derived with the newly determined 
atomic data will be lower for high-ionization objects (i.e., the Se 3+ ionic fraction will be larger 
at the expense of Se 4+ , and hence the correction for unobserved ions will be smaller) but 
quite similar for lower-excitation PNa_We contrast the Se + -Se 5+ fractional abundances from 
our current models with those of Paper I in Table [I3j In spite of the differing model input 
parameters, for most objects it can be seen that the Se 3+ ionic fraction is larger than found 


m 


Paper l typically due to a smaller Se 4+ fr action. This follows th e expe cted trend given the 


differences between the Se atomic data from Sterling & Witthoeft (20111) and those adopted 


m 


Pa per I. The primary reason that the overall Se ab undances are not drastically changed 


is that the Se 4+ fractional abundances found in Pape r ij were relatively small (no larger than 
0.32), thereby changing the Se 3+ fraction by a modest amount. 

Some of the PNe (NGC 6741, NGC 6826, NGC 6886, and NGC 7662) do not show the 
expected shift of Se 4+ from our previous models to Se 3+ in those we present here. This is 
largely due to the different model input parameters that we adopt for these objects. For 
example, we utilize a lower T e g for the central stars of NGC 6741 and NGC 6826, shifting 
the Se ionization balance to lower charge states (the larger Se 4+ fraction in NGC 6741 is due 
to a lower amount of higher-charge states than the previous model). The differing nebular 
radii of the NGC 6886 and NGC 7662 models from those of Paper l| account for the Se 
ionic fractions in those objects. The larger outer radius of the current NGC 6886 model 
shifts the Se ionization balance to lower charge states, while the smaller outer radius of the 
NGC 7662 model has the opposite effect (leading to a larger Se 4+ fraction). Finally, in our 
new IC 418 model we note that the much larger Se 2+ to Se + DR rate coefficient produces a 
larger Se + /Se 2+ ratio than found in 


Paper I. 


3.3. Remaining Atomic Data Needs 


The imperfect agreement between the modeled and observed Kr ionization balance 
demonstrates that improvements can be made to the atom i c data fo r Se and Kr . _ Ex- 


perimental PI cross sections of these ions flLu et al.l 120 06b 


Sterling et al.l l2011t lEsteves et al.l 1201 It lEsteves et al 


2012 


a; 


Lu 


2006: 


Bizau et al. 2011 


Hinojosa et al. 2012) are in¬ 


valuable for benchmarking PI cross section calculations, particularly for complex systems 
such as near-neutral n-capture elements. Comparison of the experimental and calculated PI 
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Table 13. Comparison of Se Ionic Fractions With 


PN 

Name 

Se+ 

Se 2 + 

Se 3+ 

Se 4 + 

Se 5+ 

This work 

Paper I 

This work 

Paper I 

This work 

Paper I 

This work 

Paper I 

This work 

Paper I 

IC 418 

0.34 

0.20 

0.57 

0.67 

0.09 

0.12 





IC 5117 

0.11 

0.13 

0.20 

0.20 

0.58 

0.40 

0.11 

0.26 


0.01 

NGC 6572 

0.17 

0.14 

0.27 

0.24 

0.54 

0.45 

0.02 

0.17 



NGC 6741 

0.28 

0.29 

0.19 

0.24 

0.24 

0.19 

0.26 

0.14 

0.02 

0.09 

NGC 6790 

0.08 

0.10 

0.17 

0.17 

0.64 

0.41 

0.10 

0.32 



NGC 6826 

0.08 


0.39 

0.30 

0.52 

0.61 


0.08 



NGC 6884 

0.05 

0.01 

0.24 

0.24 

0.58 

0.50 

0.12 

0.22 


0.03 

NGC 6886 

0.11 

0.01 

0.18 

0.14 

0.29 

0.31 

0.36 

0.27 

0.04 

0.18 

NGC 7027 

0.11 

0.12 

0.14 

0.17 

0.28 

0.22 

0.23 

0.15 

0.11 

0.17 

NGC 7662 

0.02 


0.13 

0.12 

0.44 

0.50 

0.40 

0.31 


0.06 


to 


Note. — Comparison of Se ionic fractions between this work and that of 


Paper I. Ionic fractions less than 1% are omitted. 
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cross sections suggest that more computationally-intensive calculations — such as R-matrix 
and/or fully relativistic treatments — would be valuable for Kr 3+ . 


Low-temperature DR cannot be accurately treated theoretically without knowledge of 
the near-threshold autoionizing state energies. Unfortunately these energies have not been 
measured experimentally for any element beyond the second row of the Periodic Table. We 
emphasize the importance of experimen tal DR rate coeffi c ient measurements conduct ed at 


facilities such as the Test Storage Ring dSavin et al. 


1999 


particularly the Cryogenic Storage Ring fiWolf et al. 


120061 ) . 


2006: Schippers et ah 201 0) and 


which can be used to measure 


DR rate coefficients for low charge-to-mass species observed in photoionized plasmas. 


Fina lly, CT rate coefficients derived through the Landau-Zener or Demkov app roxi- 
mations ( Butler &; Dalgarno 1980 : Demkovl 1964 : Swartz 1994 : Sterling fe Stand! 2011 ) are 
accurate to within a factor of three for systems with large rate coefficients, but can have 
larger errors in systems with small rate coefficients. Detail ed quantal calculations (e.g., 


Kimnra fe Lanelll989l: IZvgelman et al.lll992c IWang et al.ll2004j) can significantly improve the 


accuracy of these data, but are time- and computationally-intensive. In addition, we have 
not considered CT reactions with He, which can also affect the ionization balance. 


In a forthcoming paper, we will conduct a Monte Carlo analysis of the sensitivity of 
Se and Kr abundance determinations to atomic data uncertainties, which will highlight the 
systems and processes that most urgently require further theoretical and/or experimental 
exploration. That study will help to prioritize atomic data needs for these two elements, 
which is critical given the intensive nature of the theoretical methods and experimental 
measurements discussed. 


4. SELENIUM AND KRYPTON IONIZATION CORRECTION FACTORS 

In this section, we derive analytical ICF formulae that can be used to estimate the 
abundances of unobserved Se and Kr ions. We consider detected (or detectable) Se and Kr 
ions in the optical and K band portions of the spectrum. 

We intend for these corrections to be widely applicable, and therefore we have con¬ 
structed a large grid of Cloudy models whose properties span a wide range of central star 
temperatures and luminosities, as well as nebular densities. Specifically, we consider cen¬ 
tral star temperatures T e g- ranging from 50,000 K to 190,000 K (the range of temperatures 
of the Rauch 2003 grid), with step sizes of 5,000 K. Stellar gravities were assumed to be 
log g = 5.5 in all models. Hydrogen densities n h = 10 3 -10 5 cm -3 are considered, with step 
sizes of 3,000 cm -3 . Central star luminosities are controlled through the ionization parameter 
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U = Q(H)/ (AnR'f n nuc ), where Q(H) is the number of hydrogen-ionizing photons emitted per 
second, R in is the inner radius of the nebula, and c is the speed of light. Given the hydrogen 
density and inner radius, U is directly related to the ionizing luminosity of the central star. 
In our grids, we vary log(f/) from —3.0 to 2.0, in step sizes of 0.5. A model is computed for 
each combination of T e g, rin , and log(I/), for a total of 10,471 individual models in the grid. 

Each model is assumed to have spherical geometry and uniform density, with an inner 



we show that our derived IGF formulae are valid over a wide range of nebular compositions, 
metallicity, and assumptions regarding dust composition. 


We derive ICFs based on correlations between the Se and Kr ionic fractions and those 
of routinely detected species in the optical spectra of PNe. For each model we extracted the 
fractional abundance of all He, N, O, and Ne ions, and the first six ions of S, Cl, Ar, Se, and 
Kr. Combining the optical and K band spectra of PNe, it is possible to detect transitions 
of [Se III], [Se IV], [Kr III], [Kr IV], and [Kr V]. We searched for correlations between the 
fractional abundances of these species (and combinations thereof) and ionic fractions of He, 
N, O, Ne, S, Cl, and Ar. 


Figure [Tj depicts the strongest correlations that we found for Se and Kr ionic fractions. 
Each gray dot corresponds to the ionic fractions from a single Cloudy model. In some of the 
correlation plots, families of curves can be distinguished. These curves usually correspond 
to models with the same central star temperature, with the ionic fractions changing with 
the density and ionization parameter. 


To determine ICFs, we fit each correlation in Fig. Q]with an analytical function, using a 
X 2 minimization routine written in IDL. We considered polynomial, exponential, and loga¬ 
rithmic functional forms, attempting several types of functions (and adjusting their orders) 
to provide the best fit to each correlation. The fits are displayed as solid curves in the panels 
of Fig. m The inverses of these functions can be utilized as ICF formulae to correct for the 
abundances of unobserved Se and Kr ions in PNe. 


As in iPaper 1 we find that the Kr 2+ fractional abundance correlates strongly with the 
S 2+ and Ar 2+ ionic fractions (top left and middle panels of Fig. [[]). We derive the following 
expressions for the ICFs: 


ICF(Kr) = Kr/Kr 2+ = (-0.8774 + O.SgeSe 0 - 5735 *)" 1 , 
x = S 2+ /S; 


( 1 ) 
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0.00 


0.50 


1.00 


s 2+ / s 




0 2+ / 0 


Ar + / Ar 


0 2+ / 0 




Ar 2+ / Ar 0 2+ / 0 (Cl 2+ + Cl 3+ ) / Cl 


Fig. 1.— The strongest correlations between the fractional abundances of observable Se 
and Kr ions and those of commonly detected light element ions are plotted. These data 
are from a grid of models with solar abundances and oxygen-rich dust grains. Each gray 
dot represents the ionic fractions from one of the 10,471 models in our grid. Analytical fits 
to these correlations (see text) are shown as solid black curves, while fits from Paper I are 
shown as dashed curves if the correlation is between the same ions. Families of curves are 
seen in some of the correlations. Each of these curves generally corresponds to a single value 
of T eS , with variations in n H and U producing the changes along the curve. 
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and 

ICF(Kr) = Kr/Kr 2+ = (0.02038 + 0.4577a; + 4.752 x 10~ 7 e 13 - 96x )-\ (2) 

x = Ar 2+ /Ar. 


Equation [2] is likely to result in the most reliable correction for unobserved Kr ions, 
given the “sulfur anomaly” observed in Galactic and Magellanic Cloud PNe (Henry et al. 
2004] ; IShaw et ah 2012). A variety of causes have been sought to explain the deficit of the S 
abundance in PNe relative to H II regions and stars, but the most likely expla nati on is the 
inabi lity of current ICF schemes to accurately account for unobserved S ions flHenrv et al. 
2012; Shingles & Karakas 20131) • If this is correct, then S abundances in PNe are less accurate 
than those of Ar, which is the reason for our preference for Equation [2] 


We display the fits from Paper l| as dashed lines in Figured] Paper I found lower Kr 2+ 
fractional abundances for given S 2+ and Ar 2+ ionic fractions, meaning that their (Kr/Kr 2+ ) 
ICFs systematically overestimate Kr elemental abundances in PNe. This is in agreement 
with the hi gher Kr 2+ ionic fractions found in our models of individual PNe (§3.2) compared 


to those of 


Pa pe r I. In the case of Equation d] the disagreement with the corresponding ICF 


from [P aper I is less than 10% at the highest values of S 2+ /S (> 0.64) in our grid, but reaches 
a factor of two at the lowest values. Equation [2] shows starker disagreement with the ICF 
from Paper I: from ~ 25% for large Ar 2+ /Ar values to an order of magnitude at the lowest 
fractional abundances. 


The top right panel of Fig. Q] shows the correlation used to derive ICFs when Kr 3+ is 
the only Kr ion detected, while the middle panels show correlations for combinations of Kr 
ions. The corresponding ICF formulae are: 


ICF(Kr) = Kr/Kr 3+ = (0.06681 + 1.05a; + 0.7112a; 2 - 0.907a: 3 )- 1 , (3) 

x = Ar 3+ /Ar; 

ICF(Kr) = Kr/(Kr 2+ + Kr 3+ ) = (-2.648 + 1.047a; a8749 + 2.678e- a01552 ' T )-\ (4) 

x = 0 2+ /0; 

ICF(Kr) = Kr/(Kr 3+ + Kr 4+ ) = (0.04747 + 1.512a; - 0.5257a; 2 )- 1 , (5) 

x = Ar 3+ /Ar < 0.932; 


ICF(Kr) = Kr/(Kr 2+ + Kr 3+ + Kr 4+ ) = (0.02699 + 1.511a: - 0.5571a; 2 )- 1 , (6) 

x = 0 2+ /0. 


and 
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The upper limit to x — Ar 3+ /Ar in Equation [5] provides the range of validity for this ICF 
(above this value, the ICF is smaller than unity). None of the models in our grid achieved such 
high fractional abundances of Ar 3+ , and hence this equation should be generally applicable. 
Paper I used different reference ions than we utilize for the (Kr 2+ + Kr 3+ ) ICF, and hence a 
direct comparison is not possible. 

The bottom three panels of Fig. |T] show correlations for Se ionic fractions, including 
that of Se 2+ , whose abundance can in principle be derived from [Se III] A8854.0. The 
Se 2+ ionic abundance has not yet been reliably determined in PNe, due to the lack of 
collisional e xcitation data and because of blending with a wea k He I line at nearly the same 
wavelength (IPequignot fe Baluteaul 11994: ISharpee et al.l 120071) . From these correlations we 
find the following ICF formulae: 


and 


ICF(Se) = Se/Se 2+ = (2.411 + 0.5146a: 12 - 29 - 2.417e-°- 172;te )- 1 , 
x = Ar 2+ /Ar > 0.0344; 

ICF(Se) = Se/Se 3+ = (0.03238 + 0.3225a; + 0.4013a; 2 )- 1 , 

x = 0 2+ /0; 

ICF(Se) = Se/(Se 2+ + Se 3+ ) = (-0.3703 + O^Se 1 ’ 27 ®) -1 , 
x = (Cl 2+ + Cl 3+ )/Cl > 0.0315. 


(7) 

( 8 ) 

(9) 


The Se 3+ correlation shows the largest dispersion of these three. The situation is exacer¬ 
bated at low 0 2+ fractional abundances (< 0.25), where the Se 3+ ionic fraction has a bimodal 
distribution in our models. The bimodality of this correlation is a result of variations in the 
ionization parameter U. For small values (log (U) < —1.5), the “upper branch” tends to be 
populated, while models with larger log(C) generally populate the “lower branch.” 


Paper I also found 0 2+ to be the best tracer of the Se 3+ fractional abunda nce. Our 


models show that Se 3+ comprises a larger fractionof the Se abundance than found in Paper I. 
Equation [S] produces ICFs smaller than the Paper I values by 30-50% when 0 2+ /0> 0.5, 
and by more than an order of magnitude at the smallest 0 2+ /0 values. 


Originally, we planned to extend our grid of models to lower central star tempera¬ 
tures (down to ~30,000 K). Since the R auch (120031) grid of PN central star atmospheres 
does n ot exte n d to such low temperatures, we utilized LTE plane-parallel ATLAS models 
( Castelli fe Kurucz 2004 ) for the range T e g = 30, 000 K to 50,000 K. However, we found that 
the ionic fractions from these models do not follow the correlations depicted in Figure [D 
This suggests that either the ATLAS stellar atmosphere models do not accurately describe 
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the ionizing radiation from low-temperature PN central stars, or that our ICFs break down 
for objects with sufficiently low central star temperatures. Therefore, Equations [TH9] should 
be used with caution for PNe with cool central stars or for H II regions. However, we note 
that the Kr abundance derived with these ICFs for IC 418, the only object in our sample 
with T eS < 50,000 K (see Tabled, agrees well with the model results (see 1 J4.2H . 


4.1. Robustness of the Derived ICFs 


The grid described in the previous section was recomputed with different nebular abun¬ 
dances and dust properties, in order to test the robustness of the derived Se and Kr ICFs 
for various compositions. 

In Figure [21 we display the same correlations as shown in Figured] for a grid of mod¬ 
els that has a solar composition, but dust is not included. There is not as pronounced a 
bimodality in the Se 3+ correlation at small 0 2+ fractional abundances as seen in the grid in¬ 
cluding dust. Specifically, the “lower branch” of this correlation, at which the Se 3+ fractional 
abundance is small compared to 0 2+ , is not as highly populated in the dust-free models. 
That being said, Equation [8] falls well within the dispersion of this correlation. For other 
ions, the fits describe all of the correlations in the dust-free grid to good approximation, with 
a slightly larger dispersion in some cases. We conclude that our ICFs are robust regardless 
of whether dust is included in the models. 


The correlatio ns are s hown in Figure [3] for a grid with Cloudy’s default (C-rich) PN 
abundances flAller fc Czvzaklll983l) that does not include dust. The correlations are virtually 
identical to those found in the dust-free solar composition grid, indicating that the ICFs do 
not depend strongly on nebular composition at near-solar metallicities. Similarly, we find 
that models with Cloudy’s PN abundances that include carbon-rich dust grains (the “grain 
Orion graphite” command, where the grain properties are taken from Baldwin et al, 199111 
agree well with the correlations depicted in Figure CD 

At lower metallicities, Equations [lH9] also describe the ionic fraction correlations well. 
Figure [4] shows the fractional abundance plots of Figure [Tj for a scaled solar composition 
(without dust) in which the abundances of all elements heavier than He have been reduced 
by a factor of ten. NLTE stellar atmosphere models with halo abundances ( Rauch 2003) 
were used to describe the central star radiation in this grid. The correlations in the top 
panels, for single Kr ions, show excellent agreement with the solar metallicity grid. For 
the correlations depicted in the middle panels, the fits tend to slightly overestimate the 
fractional Kr abundances when the sum of the Kr 2+ and Kr 3+ ionic fractions is small (< 0.1- 
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Fig. 2.— The same as Figure HJ except in this grid dust is not included. The analytical fits 
to the correlations from Figure [Q are overplotted in each panel. 
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0 2+ / 0 




Ar 2+ / Ar 0 2+ / 0 (Cl 2+ + Cl 3+ ) / Cl 


Fig. 3.— The s ame as Figure CD but from a grid of models with the default Cloudy PN 
abundances flAller fe Czvzald 1 9831 ). in which C/O > 1, and no dust. The analytical fits to 
the correlations from Figure [Dare overplotted in each panel. The correlations show excellent 
agreement with those seen in the dust-free solar composition models (Figure [2]), indicating 
that the ICFs do not depend significantly on the carbon/oxygen chemistry. Similarly, results 
for C-rich nebular compositions with C-rich dust grains agree well with solar composition 
models with O-rich dust (Figure [U). 
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0.5, depending on the correlation). This will cause elemental Kr abundances derived from 
multiple ions to be slightly underestimated for some low-metallicity PNe. The Se 2+ fractional 
abundance can be underestimated by Equation [7] at low metallicity when Ar 2+ /Ar > 0.5, 
which would lead to overestimated Se abundances. 

The correlations shown in Figure [4] are very similar to those found from a grid of models 
with metal abundances reduced to 1% solar (not shown), though the discrepancies with the 
solar metallicity models described above are slightly more pronounced. Likewise, grids with 
supersolar abundances (metal abundances increased by 50% from the solar values) show 
negligible deviation from the solar abundance models with O-rich dust shown in Figure CD 

These tests show that our derived ICFs exhibit only a slight dependence on dust physics 
and metallicity. The ICFs for Kr 2+ , Kr 3+ , and Se 2+ appear to be quite robust, while those 
for combinations of Kr ions show discrepancies only in low-metallicity objects when their 
sum represents a small fraction of the total Kr abundance. We conclude that our analytical 
ICFs are valid for a wide range of nebular compositions, with errors becoming significant 
only for low-metallicity, low-ionization nebulae. 


4.2. Comparison of ICF- and Model-Derived Se and Kr Abundances 


We have computed empirical Se and Kr abundances from the optical and near-infrared 
spectra of the 15 PNe that we individually modeled (Table |TJ) , using all possible ICFs (Equa¬ 
tions [1H9]). The O, S, and Ar ionic fractions needed for the ICFs are taken from the references 
listed in Table [0 In cases where error bars were not provided for the ionic and/or elemental 
O, S, and Ar abundances, we assume a 30% uncertainty in the ionic fractions. 


Kr ionic abundances were determin ed us ing a five-level model atom calculation with 
effective collision st r engths computed bv ISchonind (119971) and transition probabilities from 
Biemont fe Hansen) (Il986a .b). Electron temperatures and densities were adopted from the 
references in Table [T] assuming uncertainties of 1000 K in T e and 20% in n e if values were 
not provided in those references. The Kr ionic abundances for each of the 15 modeled PNe 
are given in Table [14l where the error bars take into account uncertainties in line intensities 
(assumed to be 25%), electron temperature, and electron density. 

Ionic Kr abundances derived from different lines of the same emitting species agree 
within the uncertainties with the exception of Kr 2+ in IC 418 (perhaps due to aperture 
differences in the optical and near-infrared spectra, or incomplete removal of telluric OH 
contamination in A6826.7). Paper I discussed the identities of optical and near-infrared Se 
and Kr emission lines, refuting previously proposed alternative identifications. The good 
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0 2+ / 0 


Ar 3+ / Ar 0 2+ / 0 




Ar 2+ / Ar 0 2+ / 0 (Cl 2+ + Cl 3+ ) / Cl 


Fig. 4.— The same as Figure H] but with metal abundances reduced by a factor of ten from 
a solar composition. The analytical fits to the correlations from Figure |T] are overplotted 
in each panel. Small discrepancies with the fits derived from higher-metallicity grids are 
apparent for some of the correlations in the middle and bottom panels, though the Kr 2+ , 
Kr 3+ , and Se 3+ correlations appear to be robust at various metallicities. The plotted ionic 
fractions are similar to those found in a grid with metal abundances reduced by two orders 
of magnitude. 
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agreement in ionic abundances derived from different emission lines of the same Kr ion 
provides further evidence that these lines indeed arise from Kr fine-structure transitions. In 
particular, note the excellent agreement in Kr 4+ abundances derived from the A6256.1 and 
A8243.4 features in NGC 7027. These two lines arise from the same upper level, and therefore 
are independent of temperature a nd den sity. Tlnyrelative intensities validate the accuracy 
of the transition probabilities of iB iemont & Hansen (1l986an . with the small difference likely 
due to uncertainties in the line fluxes. To our knowledge, the detection of [Kr V] A8243.4 by 


Sha rpee et ah (1200711 is the first in a high-resolution PN spectrum. 


We converted the ionic abundances to elemental Kr abundances using all of the relevant 
ICFs (Table [15]). The Kr abundances found from our Cloudy models of each object are given 
in the last column of the table. Some variation is expected in the abundances derived with 
different ICFs, due to observational uncertainties (e.g., line fluxes, electron temperatures, 
and extinction coefficients). In addition, the abundances will generally be more accurate 
(i.e., the ICFs smaller) when more ions are detected. The empirical Kr abundances derived 
using the maximum number of detected Kr ions (Kr 2+ and Kr 3+ for most objects, with the 
addition of Kr 4+ in NGC 7027 and NGC 2440) agree to within ~ 0.25 dex - and often better - 
with the modeled Kr abundances. We refer to the Kr abundances derived from the maximum 
number of detected Kr ions, using Equations [4H6] as the “best empirical” Kr abundances 
since the corresponding ICFs are the smallest and least sensitive to uncertainties. Elemental 
abundances derived from a single Kr ion show more scatter (Table [T5jh as expected, though 
in many cases there is good agreement with the best empirical values. 
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Table 14. Empirical Kr Ionic Abundances in the Modeled PNe 


PN 

Name 

Kr 2+ /H+ xlO 9 
A6826 A2.199 /mi 

Kr 3 +/H+ xlO 9 
A5346 A5868 

Kr 4+ /H+ xlO 9 
A6256 A8243 

IC 418 

3.96T0.88 2.89T0.05 

0.35T0.10 0.24T0.07 


IC 2501 

0.49T0.14 

0.53T0.16 0.50T0.15 


IC 4191 

0.14T0.04 

1.80T0.62 0.98T0.29 


IC 5117 

1.35T0.34 

5.92T1.50 4.90T1.19 


NGC 2440 


0.91T0.20 0.88T0.23 

0.25T0.05 

NGC 5315 

0.90T0.26 

1.19T0.33 


NGC 6369 

1.39T0.38 

4.15T1.24 4.13T1.21 


NGC 6572 

1.45T0.55 0.88T0.21 

2.89T0.93 2.74T0.97 


NGC 6741 

2.03T0.35 1.57T0.40 

2.60T0.74 1.73T0.56 


NGC 6790 

0.64T0.09 

1.99T0.56 1.36T0.32 


NGC 6826 

1.77T0.43 

2.05T0.75 2.77T0.69 


NGC 6884 

1.12T0.24 

5.22T1.62 5.07T0.96 


NGC 6886 

1.80T0.37 

3.45T1.01 2.77T0.79 


NGC 7027 

2.55T0.67 2.36T0.41 

5.68T1.28 5.43T1.21 

1.07T0.29 0.98T0.27 

NGC 7662 

0.63T0.13 

5.16T1.42 4.63T0.63 



Note. — Comparison of Kr ionic abundances derived from optical and near-infrared 
emission lines. See Table |T] for references to the spectroscopic data utilized. 




















Table 15. Comparison of ICF- and Model-Derived Kr Abundances 


PN 

12+ Log(Kr/H) 

12+ Log(Kr/H) 

12+ Log(Kr/H) 

12+ Log(Kr/H) 

12+ Log(Kr/H) 

12+ Log(Kr/H) 

12+ Log(Kr/H) 

Name 

ICF Eqn[2] 

ICF Eqn[T| 

ICF Eqn|3] 

ICF Eqng] 

ICF Eqn[5] 

ICF Eqn[6] 

Model 

IC 418 

4.20+0.15 

3.87+0.20 

a 

4.06+0.18 



3.93+0.23 

IC 2501 

3.01+0.33 

3.02+0.24 

3.55+0.12 

3.02+0.09 



3.12+0.13 

IC 4191 

2.83+0.19 

2.81+0.23 

3.53+0.17 

3.20+0.09 



3.22+0.30 

IC 5117 

3.64+0.21 

3.92+0.22 

4.07+0.15 

3.90+0.12 



3.74+0.10 

NGC 2440 



3.44+0.18 


3.44+0.11 


3.44+0.10 

NGC 5315 

3.37+0.30 

3.24+0.13 

3.68+0.17 

3.36+0.10 



3.42+0.28 

NGC 6369 

3.54+0.31 

3.73+0.23 

4.14+0.15 

3.77+0.10 



3.94+0.18 

NGC 6572 

3.48+0.29 

3.66+0.22 

3.99+0.16 

3.60+0.08 



3.51+0.15 

NGC 6741 

3.78+0.21 

3.81+0.23 

3.81+0.16 

3.95+0.19 



3.69+0.20 

NGC 6790 

3.35+0.16 

3.42+0.19 

3.50+0.15 

3.37+0.08 



3.34+0.10 

NGC 6826 

3.53+0.34 

3.61+0.23 

4.16+0.13 

3.64+0.09 



3.63+0.11 

NGC 6884 

3.61+0.17 

3.86+0.20 

3.98+0.16 

3.82+0.10 



3.69+0.15 

NGC 6886 

3.82+0.20 

3.94+0.16 

3.94+0.16 

3.82+0.20 



3.70+0.20 

NGC 7027 

3.94+0.16 

4.10+0.20 

4.02+0.15 

3.97+0.20 

4.09+0.11 

4.01+0.10 

4.14+0.11 

NGC 7662 

3.49+0.16 

4.03+0.21 

3.87+0.15 

3.96+0.19 



3.74+0.11 

PN 

ICF 

ICF 

ICF 

ICF 

ICF 

ICF 


Name 

Eqn[2] 

Eqn[T] 

Eqn[3] 

Eqn(4] 

Eqn[5] 

Eqn[6] 


IC 418 

4.60+1.36 

2.14+0.18 


3.12+1.22 




IC 2501 

2.10+1.21 

3.02+0.24 

6.97+1.19 

1.04+0.22 




IC 4191 

4.80+1.41 

4.63+1.79 

2.44+0.69 

1.05+0.22 




IC 5117 

3.21+1.11 

6.21+2.37 

2.18+0.62 

1.17+0.30 




NGC 2440 



3.10+0.83 


2.82+0.72 



NGC 5315 

2.58+1.36 

1.92+0.07 

4.01+0.98 

1.09+0.25 




NGC 6369 

2.48+1.36 

3.83+1.50 

3.36+0.88 

1.08+0.24 




NGC 6572 

2.58+1.36 

3.93+1.53 

3.50+0.90 

1.00+0.19 




NGC 6741 

3.37+1.11 

3.56+1.40 

2.99+0.81 

2.26+0.92 




NGC 6790 

3.54+1.13 

4.14+1.61 

1.88+0.54 

1.00+0.19 




NGC 6826 

1.90+1.15 

2.29+0.94 

6.01+1.16 

1.04+0.22 




NGC 6884 

3.66+1.14 

6.45+2.46 

1.84+0.52 

1.07+0.24 




NGC 6886 

3.60+1.13 

5.23+2.01 

2.78+0.76 

1.34+0.57 




NGC 7027 

4.00+1.22 

5.78+2.21 

2.24+0.64 

1.35+0.58 

2.09+0.54 

1.27+0.27 


NGC 7662 

4.91+1.44 

16.9+6.8 

1.52+0.41 

1.64+0.69 
























Note. — Comparison of Kr elemental abundanc es (top) derived with different ICFs and from Cloudy models. The solar Kr abundance is 
12 + Log(Kr/H) = 3.26 =L 0.06 dAsplund et alj|2009h . The values of the ICFs are given in the lower portion of the table. 


a Eq uation (3) is not used to determine the Kr abundance of IC 418, since Ar 3 + was not detected in 
d2004h . 


the optical spectrum used by lPottasch et ahl 
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In Figure [5l we plot the Kr abundances determined from single Kr ions (using Equa¬ 
tions |T] and [2] for Kr 2+ , and Equation [3] for Kr 3+ ) against the best empirical Kr abundance. 
Abundances derived with Equation [Tj agree very well with the best empirical values. This 
i s sur prising in light _of the uncertainty regarding total elemental S abundances in PNe 
( Henry et al.1 20041 . 2012), which can affect the value of the ICF. Equation [2] underestimates 
the gaseous Kr abundance in some of the PNe, though the agreement with the best empir¬ 
ical value is typically quite good (less than 0.2 dex). We nevertheless regard Equation [2] to 
be the most reliable ICF when only Kr 2+ emission is detected, both because of the sulfur 
abundance anomaly in PNe and because the S 2+ abundance (needed for Equation [T]) is often 
derived from the strongly temperature-sensitive A6312.1 line. 


There is a tendency for Equation E] to produce larger Kr abundances than the best 
empirical estimate. This is not a systematic trend, as in some cases the elemental abundances 
are in good agreement (e.g., NGC 6741, NGC 6790, NGC 6886, NGC 7027, NGC 7662). 
Equation [3] is largest (and hence most sensitive to uncertainties) when the Ar 3+ /Ar fraction 
is small. Indeed, many of the larger disparities are seen for objects in which Ar 3+ /Ar < 0.20 
(IC 2501, NGC 5315, NGC 6369, NGC 6572, NGC 6826). Fortunately, these objects are 
low- or moderate-excitation PNe (with T e g typically between 50,000 K and 60,000 K) in 
which Kr 2+ can be readily detected. Therefore, the ICF from Equation [3] should be used 
with caution in such objects, in which Equations [2] and |4] are preferable if [Kr III] emission 
is detected. 


The only other object with a discrepant Kr abun dance from Equation's IC 4191. This 
is most likely due to observational error, as ISharpee et ah (120071) flag [Kr III] A6826.7 as a 
marginal detection. For this object, the disagreement in Kr abundances from the different 
ICFs is likely due to the uncertain Kr 2+ abundance. 

Observational errors in the O, S, and Ar ionic fractions also affect the Kr abundances. 
For example, the NGC 7662 Kr abundances derived from Equations [2] and [3] differ by 0.38 dex 
when using ionic Ar fractions from [Liu et ah (]2004af) . but only by 0.15 dex (3.68 and 3.83, 
respectively) if those from iHvung fe Allen ( 19971) are used. In addition, various ICF schema 


exist for computing Ar abun dances (e.g.. lKingsburgh fe Barlowlll994J:lKwitter &; Henrvll2001 


Delgado-Inglada et al.l l2014h which can also lead to differences in Ar ionic fractions. These 
effects also contribute to the scatter in Kr abundances derived with the various ICFs. 

Se 3+ /H + abundances were derived from the 2.287 /jrn line, using data from the near- 
infrared references in Table [0 Since the ground configuration of Se 3+ consists of a 2 P term, 
a two-level model atom is sufficient for deriving the ionic abundance (the first excited term 
has an energy ~10 eV higher tha n the ground configuration). We utilized the transition 
probability of Biemont & Hansen (jl987l ) and the collision strength calculated by K. Butler 
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Fig. 5.— A comparison of the Kr abundances derived from Kr 2+ lines (Equations [T] and [2] 
filled and open triangles respectively) and from Kr 3+ lines (Equation [3] x symbols) with the 
best empirical Kr abundances (derived from all detected Kr ions). The dashed line denotes 
perfect agreement, and is not a fit to the data. 
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(2007, private communication). 

Table [TH] displays the calculated Se 3+ abundances, as well as the empirical and modeled 
elemental abundances. The ICF- and model-derived abundances agree to within the uncer¬ 
tainties for all but NGC 6886, and are within 0.15 dex of each other with the exceptions of 
NGC 6886 and NGC 7027. The level of agreement seen is consistent with observational and 
modeling uncertainties, and with the dispersion in the ionic fraction correlations of Fig. [0 


This analysis shows that our model-derived ICF formulae produce Se and Kr elemental 
abundances that are consistent within the uncertainties with abundances derived from mul¬ 
tiple Kr ions and with Cloudy models of the individual PNe. We now apply our ICFs to the 
K band PN survey of Paper II. 


5. APPLICATION OF ICFS TO PLANETARY NEBULA SE AND KR 

ABUNDANCES 

In this section, we utilize Equations [H [2j and [8] to re-evaluate Se and Kr enrichments in 
PNe. We consider the K band spectra of 120 PNe from Paper II, and examine whether the 
correlations found in that study are affected by the changes to the Se and Kr abundances 
due to using updated atomic data. 


We aim to make the comparison with Paper II as direct as possible, and therefore utilize 
the same ionic abundances, electron temperatures, and densities as before. We change only 
the Se and Kr ICFs in or der to isolate the effects of the upd ated ato mic data. The solar 


abundances of lAsplund et al.l (120051) are adopted, fo 


lations of the solar composition (e.g., Asplund et al¬ 


lowin g iPaper III . More recent compi- 


2009) adopt a si gnific antly larger Ar 


abund ance, in better agreement with the studies of lLoddersI (120031.120081) than lAsplund et al. 
(2005). We discuss how the revised solar abundances would affect the Se and Kr enrichment 
factors in the following text and tables. 

We present the Se and Kr abundances relative to the reference elements O and Ar, in 
order to suppress effects due to intrinsic dispersion in metallicity as well as systematic effects 
such as abunda nce gradients with radial Galactocentri c distance or vertical distance fro m the 
Galact ic plane ( Maciel et al. 2006 : Henry et al. 2010 : Stanghellini fe Havwood 2010 b Fol¬ 
lowing IPaper I], we determine Se and Kr abundances relative to Ar in Type I PNe and to O 
in other objects. Type I PNe are believed to be desce ndants of more massive progenitor stars 
(> 3- 4 M 0 ), based on their He and N enrichments (|Peimbertlll978l: iKingsburoh fc Harlow 
1994b small Galactic scale height, high nebular and stellar masses, and peculiar velocities 
(ITorres-Peimbert &; PeimbertlIl997l: ICorradi &; SchwarzlIl995l: iGornv et al.lll997l: IPena et al. 
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Table 16. Comparison of ICF- and Model-Derived Se Abundances 


PN 

Name 

n(Se 3+ )/ 

n(H+)xl0 9 

ICF(Se) 

12+ Log(Se/H) 

ICF 

12+ Log(Se/H) 
Model 

IC 418 





IC 2501 





IC 4191 





IC 5117 

2.15+0.23 

1.83+0.81 

3.59+0.22 

3.52+0.05 

NGC 2440 





NGC 5315 





NGC 6369 

2.05+0.42 

1.60+0.64 

3.51+0.21 

3.66+0.14 

NGC 6572 

0.68+0.04 

1.40+0.50 

2.98+0.16 

3.01+0.03 

NGC 6741 

0.58+0.14 

4.99T1.75 

3.46+0.20 

3.31+0.11 

NGC 6790 

0.55+0.06 

1.38T0.48 

2.88+0.17 

2.91+0.05 

NGC 6826 

0.47+0.14 

1.50T0.57 

2.85+0.23 

2.95+0.13 

NGC 6884 

1.50+0.22 

1.57T0.62 

3.37+0.20 

3.40+0.06 

NGC 6886 

0.99+0.11 

2.29T1.01 

3.35+0.21 

3.70+0.08 

NGC 7027 

1.26+0.08 

2.33T1.03 

3.47+0.21 

3.68+0.08 

NGC 7662 

1.44+0.32 

3.16T1.29 

3.66+0.22 

3.53+0.10 


Note. — Comparison of Se abundances derived from ICFs and from 
Cloudy models. The solar Se abundance is 12 + Log(Se/H) = 3.34 ± 0.03 
( Asplund et al. 2009 b 
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2013 ). The primary reason for using different metallicity indicators is that the_^Ar/0 ratio 


is a factor of two larger in Type I PNe than in non-Type I objects of Paper IPs sample 


Similar trends were found for S/O and Cl/O, suggesting that uncertainties in the Ar abun¬ 


dances are not to blame. In Paper II], we speculated that this disparity in Ar/O ratios for 


the two populations could be due to ON-burning at the base of the convective envelope 

(i.e., hot bottom burning or HBB; 


Ventura & D’Antona 

2005a 

b; 

Karakas et al. 

2006) 


nucleosynthesis showed that the overall O abundance decreases only modestly (by ~ 0.06- 
0.15 dex, depending on the assumed s olar abundances) as a result of HBB for solar-met allicity 


stars with initial masses 6-6.5 M 0 ( Karakas et al. 2009 : Karakas fc Lattanzio 2014h. Fur¬ 


thermore, S, Cl, and Ar are not significantly affecte d by AGB nucleosynthesis (IKarakas et al. 


2009 : Henry et al. 20121: Shingles fe Karakas 2013 b Therefore the cause of the difference in 
the Ar/O (and S/O and Cl/O) ratios for Type I and non-Type I PNe is unresolved at this 
time. Whether it is due to systematic uncertainties in the abundance determinations of 
Type I PNe or to an unknown nucleosynthetic or mixing process in intermediate-mass AGB 
stars, we adopt Ar as the reference element in Type I PNe. Oxygen abundances are generally 
more accurate in PNe, and hence we use O as the metallicity reference in non-Type I PNe; 
the results below are not significantly affected if Ar is used as the reference element for the 
full sample. 

The Se and Kr ionic abundances used in our abundance determinations were derived 


previously and can be found in Table 9 of Paper II. Likewise, the electron temperatures, 
densities, and ionic and elemental O, S, and Ar abundances for the sample are given in their 
Tables 7 and 8. We use our new Equation [2] to derive Kr abundances, except when Ar ionic 
or elemental abundances are not available for a PN, or if Equation [2] gives an anomalously 
large ICF compared to Equation [0 Se and Kr abundances derived with the new ICFs are 
given in Table [171 relative to H, O, and Ar. The [Kr/O] and [Se/O] values for Type I PNe 
are enclosed in parentheses to indicate that they may not be reliable. If the solar abundances 
of Asplund et ah (.2009) are adopted, the listed [Kr/O] and [Kr/Ar] values would increase 
by 0.05 and 0.24 dex, respectively, while [Se/O] and [Se/Ar] would increase by 0.02 and 
0.21 dex. 
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Table 17. Se and Kr Abundances Relative to H, 0, and Ar 


Object 


Name 

[Kr/H] 

[Kr/0] a ’ 6 

[Kr/Ar] b 

[Se/H] 

[Se/0] a ’ b 

[Se/Ar] b 

Ref.® 

BD+30°3639 

<0.96 f 

< 0.96 

< 0.43 

<0.55 

< 0.54 

< 0.01 

IR2, 7 


<0.89 f 

< 0.89 

< 0.35 




IR4, 7 


0.22=t0.21 f 

0.22±0.23 

-0.324=0.23 




IR6, 7 

Cn 3-1 

<0.73 

< 0.76 

< 0.82 

<0.71 

< 0.74 

< 0.80 

IR1, 46 

DdDm 1 

<1.06 

< 1.67 

< 2.08 

<0.24 

< 0.85 

< 1.26 

IR1, 46 

Hb 4 

<0.58 

(< 0.56) 

< 0.26 

0.144=0.27 

(0.124=0.30) 

-0.184=0.34 

IR1, 2 

Hb 5 

0.364=0.18 

(0.324=0.20) 

-0.244=0.20 

-0.484=0.24 e 

(-0.524=0.26)® 

-1.084=0.26® 

IR1, 37 

Hb 6 

<0.70 

(< 0.65) 

< 0.18 

0.124=0.34 

(0.074=0.41) 

-0.404=0.41 

IR1, 2 

Hb 7 




<0.42 

< 0.43 

< 0.58 

IR1, 41 

Hb 12 

0.394=0.15 

0.71±0.18 

0.27±0.18 

-0.384=0.18 

-0.064=0.21 

-0.504=0.21 

IR5, 18 

He 2-459 

<1.25 f 

(< 2.06) 


<0.68 

(< 1.49) 


IR1, 14 

Hu 1-1 

<1.12 f 

< 1.22 

< 2.10 

<0.36 

< 0.46 

< 1.34 

IR1, 46 

Hu 1-2 

<0.87 

(< 1.45) 

< 1.09 

<0.06 

(< 0.64) 

< 0.27 

IR1, 29 

Hu 2-1 

0.514=0.30 

0.66±0.32 

0.91±0.32 

<-0.71 

< -0.56 

< -0.31 

IR1, 46 

IC 351 

<1.19 

< 1.45 

< 1.42 

0.104=0.28® 

0.364=0.31® 

0.334=0.31® 

IR1, 46 

IC 418 

0.50=t0.20 f 

0.62±0.22 

0.43±0.22 

<-0.32 

< -0.21 

< -0.40 

IR1, 39 


0.44±0.46 e ^ 

0.55±0.49 e 

0.36±0.49 e 




IR3, 39 

IC 1747 

<1.06 

< 1.14 

< 1.16 

0.374=0.22 

0.454=0.24 

0.464=0.24 

IR1, 46 

IC 2003 

<0.95 

< 1.17 

< 1.18 

0.334=0.21 

0.554=0.24 

0.564=0.24 

IR1, 46 


<0.95 

< 1.17 

< 1.18 

0.454=0.22 

0.674=0.24 

0.684=0.24 

IR2, 46 





0.394=0.31® 

0.614=0.34® 

0.624=0.34® 

IR3, 46 

IC 2149 

<-0.06 

< 0.14 

< 0.22 

<-0.51 

< -0.30 

< -0.23 

IR1, 13 

IC 2165 

<0.50 

< 0.73 

< 0.60 

-0.104=0.23 

0.134=0.25 

-0.004=0.25 

IR1, 39 


<0.59 

< 0.82 

< 0.69 

0.044=0.24 

0.274=0.26 

0.144=0.26 

IR2, 39 

IC 3568 

<0.28 

< 0.55 

< 0.57 

-0.564=0.22® 

-0.294=0.24® 

-0.274=0.24® 

IR1, 29 

IC 4593 

<0.51 

< 0.48 

< 0.43 

<-0.06 

< -0.10 

< -0.14 

IR1, 16 

IC 4634 

<0.17 

< 0.26 

< 0.35 

-0.414=0.18 

-0.324=0.21 

-0.234=0.21 

IR1, 21 

IC 4732 

<0.82 

< 1.06 

< 1.02 

0.024=0.39® 

0.264=0.47® 

0.224=0.47® 

IR1, 2 

IC 4846 

<0.56 

< 0.72 

< 0.78 

<-0.38 

< -0.22 

< -0.16 

IR1, 46 

IC 49976 

<-0.11 

< -0.53 

< 0.46 

<-0.67 

< -1.09 

< -0.10 

IR1, 20 


<-0.22 

< -0.64 

< 0.36 

<-0.90 

< -1.32 

< -0.32 

IR2, 20 

IC 5117 

0.384=0.16 

0.45±0.19 

0.35±0.19 

0.264=0.15 

0.334=0.18 

0.234=0.18 

IR7, 22 


0.484=0.28® 

0.55±0.30 e 

0.45±0.30® 

0.324=0.16 

0.394=0.19 

0.304=0.19 

IR2, 22 


0.604=0.14 

0.67±0.17 

0.574=0.17 




IR4, 22 

IC 5217 

<0.50 

< 0.66 

< 0.58 

-0.094=0.20 

0.074=0.23 

-0.014=0.23 

IR1, 46 

J 320 

<0.78 

< 1.04 

< 0.98 

<-0.11 

< 0.16 

< 0.10 

IR1, 31 

J 900 

0.514=0.30® 

0.62±0.31 e 

0.61±0.32 e 

0.324=0.25 

0.434=0.26 

0.424=0.27 

IR1, 28 

K 3-17 

>0.25 

> 0.05 


0.434=0.23 

0.244=0.28 


IR1, 25 

K 3-55 

>0.20 

> -0.04 


0.244=0.28 

0.014=0.34 


IR1, 24 

K 3-60 

1.164=0.39 

1.19±0.75 

0.924=0.75 

0.584=0.60 

0.614=0.80 

0.344=0.80 

IR1, 2 


1.134=0.39 

1.16±0.72 

0.894=0.72 

0.584=0.64 

0.614=0.84 

0.344=0.84 

IR6, 2 

K 3-61 

<1.42 

(< 1.39) 

< 1.08 

0.364=0.34 

(0.334=0.36) 

0.024=0.36 

IR1, 2 

K 3-62 

>0.04 



>-0.36 



IR6 

K 3-67 




-0.354=0.33 

0.054=0.40 

-0.174=0.40 

IR6, 2 

K 4-48 

<0.95 

< 0.94 

< 1.40 

<0.76 

< 0.75 

< 1.20 

IR1, 9 
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Table 17—Continued 


Object 

Name 

[Kr/H] 

[Kr/0] a ’ 6 

[Kr/Ar] b 

[Se/H] 

[Se/0] a ’ i> 

[Se/Ar] b 

Ref.® 


0.734=0.34 

0.714=0.35 

1.174=0.38 




IR6, 9 

M 1-1 

<1.84 

< 1.80 

< 1.69 

<1.09 

< 1.05 

< 0.95 

IR1, 3 

M 1-4 

<0.94 

< 1.10 

< 0.94 

0.304=0.34 

0.464=0.41 

0.304=0.41 

IR1, 2 


0.754=0.21 

0.914=0.28 

0.754=0.28 

0.314=0.37 

0.474=0.45 

0.314=0.45 

IR6, 2 

M 1-5 

0.414=0.38 

0.964=0.44 

0.594=0.44 

<-0.78 

< -0.23 

< -0.60 

IR1, 31 

M 1-6 

0.12=t0.24 f 

0.434=0.29 

-0.324=0.29 

<0.09 

< 0.41 

< -0.34 

IR1, 8 


0.05=t0.24 f 

0.374=0.29 

-0.384=0.29 




IR6, 8 

M 1-9 

<0.94 

< 1.28 

< 1.12 

<-0.06 

< 0.28 

< 0.12 

IR1, 8 

M 1-11 

0.794=0.34 

1.714=0.38 

1.864=0.39 

0.404=0.26® 

1.324=0.70® 

1.474=0.75® 

IR1, 9 


0.794=0.35 

1.714=0.39 

1.854=0.39 

-0.094=0.37 

0.834=0.41 

0.974=0.41 

IR6, 9 

M 1-12 

0.464=0.37 

0.684=0.38 

1.324=0.41 

<0.63 

< 0.85 

< 1.49 

IR1, 27 


0.414=0.37 

0.634=0.38 

1.274=0.41 




IR6, 27 

M 1-14 

<0.18 f 

< 0.54 

< 0.05 

<0.04 

< 0.40 

< -0.09 

IR1, 8 

M 1-16 

<0.63 

(< 0.91) 

< 0.55 

<-0.08 

(< 0.20) 

< -0.16 

IR1, 32 

M 1-17 

0.944=0.39 

0.744=0.45 

0.604=0.45 

0.254=0.30 

0.054=0.35 

-0.094=0.35 

IR1, 8 

M 1-20 




-0.104=0.34 

-0.064=0.40 

0.034=0.40 

IR6, 40 

M 1-25 

0.634=0.40 

0.544=0.46 

0.154=0.46 

<-0.02 

< -0.11 

< -0.50 

IR1, 14 

M 1-31 

<0.57 

< 0.74 d 

< 0.33 d 

0.044=0.44® 

0.214=0.51 d,e 

-0.204=0.64 d,e 

IR1, 14 

M 1-32 

1.17=t0.34 f 

1.174=0.35 

0.464=0.35 

0.514=0.23 

0.514=0.28 

-0.204=0.28 

IR1, 14 

M 1-35 

<1.09 

(< 1.27) 

< 0.64 

0.144=0.48 

(0.324=0.59) 

-0.314=0.59 

IR1, 2 

M 1-40 

0.224=0.24 

(0.394=0.24) 

-0.204=0.29 

-0.024=0.08® 

(0.154=0.08)® 

-0.454=0.16® 

IR1, 15 

M 1-46 

<1.08 

< 0.87 

< 0.73 

<1.17 

< 0.96 

< 0.82 

IR1, 14 

M 1-50 

<0.86 

< 0.70 

< 0.67 

-0.044=0.40® 

-0.214=0.46® 

-0.234=0.46® 

IR1, 31 

M 1-51 

1.024=0.36 

0.764=0.37 

0.664=0.40 

0.534=0.36 

0.274=0.42 

0.164=0.56 

IR1, 14 

M 1-54 

<1.05 f 

(< 0.94) 


<0.44 

(< 0.33) 


IR1, 31 

M 1-57 

0.544=0.28 

(0.394=0.33) 

0.054=0.33 

-0.044=0.37 

(-0.194=0.43) 

-0.534=0.43 

IR1, 31 

M 1-58 

<1.03 

< 0.93 

< 0.74 

0.344=0.30 

0.244=0.35 

0.054=0.83 

IR1, 9 

M 1-60 

<0.74 

(< 0.62) 

< 0.24 

0.244=0.51 

(0.124=0.60) 

-0.254=0.59 

IR1, 14 

M 1-61 

<0.31 

< 0.48 

< 0.32 

-0.214=0.39 

-0.044=0.45 

-0.204=0.45 

IR1, 14 

M 1-71 

0.564=0.26 

0.464=0.31 

0.244=0.31 

-0.094=0.43 

-0.194=0.49 

-0.414=0.49 

IR1, 14 

M 1-72 

<0.17 

< 0.90 

< 0.78 

<0.33 

< 1.06 

< 0.94 

IR1, 27 

M 1-74 

<0.38 

< 0.44 

< 0.32 

-0.644=0.44® 

-0.584=0.46® 

-0.704=0.46® 

IR1, 46 





-0.414=0.23 

-0.354=0.25 

-0.474=0.25 

IR6, 46 

M 1-75 

<1.41 f 

(< 1.28) 

< 0.75 

<0.93 

(< 0.80) 

< 0.27 

IR1, 2 

M 1-80 

<1.09 f 

< 1.16 

< 1.14 

0.454=0.29 

0.524=0.37 

0.514=0.36 

IR1, 2 

M 2-2 

<1.30 

< 1.53 

< 1.53 

<0.22 

< 0.45 

< 0.45 

IR1, 2 

M 2-31 

<0.75 

< 0.73 

< 0.56 

0.054=0.43 

0.034=0.50 

-0.144=0.50 

IR1, 40 

M 2-43 

0.814=0.35 

0.844=0.39 

0.824=0.39 

0.064=0.32 

0.094=0.37 

0.084=0.37 

IR1, 14 

M 2-48 

<1.42 

(< 1.53) 

< 0.73 

<0.33 

(< 0.44) 

< -0.36 

IR1, 30 

M 3-15 

<0.50 

< 0.80 

< 0.54 

-0.294=0.38 

0.014=0.44 

-0.254=0.58 

IR1, 14 

M 3-25 

<0.56 

(< 0.39) 

< 0.22 

0.394=0.39 

(0.224=0.45) 

0.054=0.45 

IR1, 9 

M 3-28 

<2.10 f 

(< 2.03) 


<0.42 

(< 0.35) 


IR1, 25 

M 3-35 




>-0.75 



IR1, 5 


M 3-41 


(< 2 . 10 ) 


<1.50 (< 2.07) 


IR1, 40 


M 4-18 


IR1, 10 


<1.53 


< 1.72 


< 1.69 






Table 17—Continued 


Object 


Name 

[Kr/H] 

[Kr/0] a ’ b 

[Kr/Ar] b 

[Se/H] 

[Se/0] a ’ 6 

[Se/Ar] b 

Ref.® 

Me 1-1 

<0.81 

(< 0.91) 

< 0.63 

<-0.07 

(< 0.03) 

< -0.25 

IR1, 42 

Me 2-1 




0.074=0.21 

0.014=0.23 

0.054=0.23 

IR6, 43 

Me 2-2 

<0.43 

(< 0.76) 

< 0.84 

<-0.42 

(< -0.09) 

< -0.01 

IR1, 46 

NGC 40 

1.09±0.19 f 

1.03±0.21 

0.74=1=0.21 




IR4, 38 


1.13±0.20 f 

1.07±0.23 

0.78=1=0.23 




IR3, 38 

NGC 1501 

<1.76 

(< 1.89) 

< 1.83 

<0.95 

(< 1.07) 

< 1.02 

IR1, 12 

NGC 2392 

<0.92 

< 1.00 

< 0.88 

<0.17 

< 0.26 

< 0.14 

IR1, 16 

NGC 2440 

<0.52 

(< 0.61) 

< 0.20 

<-0.30 

(< -0.22) 

< -0.62 

IR2, 6 

NGC 3242 

<0.73 

< 0.87 

< 0.92 

0.014=0.21 

0.154=0.23 

0.204=0.23 

IR1, 44 


<0.16 

< 0.30 

< 0.34 

0.014=0.19 

0.154=0.21 

0.204=0.21 

IR2, 44 





0.054=0.23 

0.194=0.25 

0.244=0.25 

IR3, 44 





-0.144=0.25 

-0.004=0.27 

0.044=0.27 

IR3, 44 

NGC 6210 

<0.18 

< 0.20 

< 0.24 

-0.214=0.18 

-0.204=0.21 

-0.154=0.21 

IR2, 29 





-0.274=0.24 

-0.254=0.27 

-0.214=0.27 

IR3, 29 

NGC 6302 

<0.66 

(< 0.95) 

< 0.06 

<-0.14 

(< 0.16) 

< -0.74 

IR2, 33 

NGC 6369 

<1.25 f 

< 1.45 

< 0.81 

0.434=0.34 

0.634=0.42 

-0.014=0.42 

IR1, 2 

NGC 6439 

<1.20 f 

« 1-21) 

< 0.93 

0.254=0.27 

(0.264=0.30) 

-0.024=0.30 

IR1, 2 

NGC 6445 

1.09±0.28 e 

0.704=0.33® 

0.77±0.33 e 

1.024=0.26 

0.634=0.31 

0.704=0.31 

IR1, 27 

NGC 6537 

<0.25 f 

(< 0.64) 

< -0.19 

-0.314=0.20 

(0.094=0.22) 

-0.744=0.22 

IR1, 34 


<0.45 f 

(< 0.85) 

< 0.02 

<-0.47 

(< -0.08) 

< -0.90 

IR2, 34 

NGC 6543 




-0.164=0.28 

-0.364=0.30 

-0.494=0.30 

IR3, 45 

NGC 6567 

<0.34 

< 0.70 

< 0.85 

-0.084=0.18 

0.284=0.21 

0.434=0.21 

IR1, 19 

NGC 6572 

0.08±0.28 

0.184=0.30 

0.07=1=0.30 

-0.354=0.14 

-0.244=0.17 

-0.364=0.17 

IR1, 29 


0.19±0.31 e 

0.294=0.33® 

0.18±0.33 e 

-0.394=0.16 

-0.294=0.18 

-0.404=0.18 

IR2, 29 


0.15±0.49 e 

0.254=0.53® 

0.144=0.53® 

-0.394=0.18 

-0.294=0.21 

-0.404=0.21 

IR3, 29 


0.31±0.39 e 

0.414=0.41® 

0.30±0.41 e 

-0.284=0.16 

-0.184=0.19 

-0.304=0.19 

IR3, 29 

NGC 6578 

<1.03 

< 0.94 

< 0.71 

0.314=0.28 

0.224=0.92 

-0.014=0.92 

IR1, 2 

NGC 6629 

0.96±0.26 

1.02=1:0.33 

0.544=0.33 

-0.324=0.68® 

-0.264=1.00® 

-0.744=1.00® 

IR1, 2 

NGC 6644 

0.45±0.26 e 

0.544=0.28® 

0.514=0.28 e 

-0.134=0.21 

-0.044=0.24 

-0.074=0.24 

IR1, 4 

NGC 6720 

1.14±0.31 e 

1.004=0.32® 

0.784=0.32 e 




IR3, 29 

NGC 6741 

0.44±0.19 

0.264=0.22 

0.254=0.22 

0.144=0.19 

-0.054=0.22 

-0.064=0.22 

IR1, 29 

NGC 6751 

<1.65 

< 1.53 

< 1.30 

<1.07 

< 0.94 

< 0.71 

IR1, 26 

NGC 6778 

<0.85 

(< 1.02) 

< 0.60 

<0.34 

(< 0.51) 

< 0.08 

IR1, 1 

NGC 6790 

<-0.05 

< 0.21 

< 0.37 

-0.444=0.14 

-0.184=0.17 

-0.024=0.17 

IR1, 29 

NGC 6803 

<0.23 

< 0.19 

< 0.03 

-0.294=0.18 

-0.334=0.20 

-0.504=0.20 

IR1, 46 


-0.02±0.43 e 

-0.064=0.44® 

-0.234=0.44® 

-0.224=0.20 

-0.264=0.22 

-0.434=0.22 

IR3, 46 

NGC 6804 




<0.82 

< 0.98 


IR1, 2 

NGC 6807 

<0.65 f 

< 0.75 

< 1.35 

<-0.55 

< -0.46 

< 0.14 

IR1, 46 

NGC 6818 

<0.50 

< 0.47 

< 0.24 

0.344=0.20 

0.314=0.22 

0.084=0.22 

IR1, 35 

NGC 6826 

<0.37 

< 0.51 

< 0.56 

-0.494=0.22 

-0.354=0.25 

-0.304=0.25 

IR1, 29 

NGC 6833 

<0.26 

< 0.79 

< 0.63 

<-0.53 

< -0.00 

< -0.17 

IR1, 46 

NGC 6879 

<0.57 

< 0.71 

< 0.52 

0.014=0.18 

0.154=0.21 

-0.044=0.21 

IR1, 46 

NGC 6881 

0.74±0.23 e 

0.664=0.30 d > e 

0.274=0.30 d,e 

0.304=0.25 

0.224=0.32 d 

-0.174=0.32 d 

IR1, 2 

NGC 6884 

<0.43 

< 0.48 

< 0.35 

0.024=0.18 

0.074=0.21 

-0.064=0.21 

IR1, 29 

NGC 6886 

0.50±0.17 

0.35=1=0.19 

0.364=0.19 

0.034=0.20 

-0.124=0.22 

-0.114=0.22 

IR1, 36 
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Table 17—Continued 


Object 

Name 

[Kr/H] 

[Kr/0] a ’ b 

[Kr/Ar] b 

[Se/H] 

[Se/0] a ’ 6 

[Se/Ar] b 

Ref. c 

NGC 6891 

<0.54 

< 0.59 

< 0.57 

<-0.29 

< -0.23 

< -0.25 

IR1, 46 

NGC 6905 

<1.54 

< 1.33 

< 1.62 

<1.22 

< 1.01 

< 1.30 

IR1, 26 

NGC 7009 




-0.09±0.32 e 

-0.01±0.34 e 

-0.33±0.34 e 

IR3, 17 

NGC 7026 

<0.23 

< 0.15 

< 0.02 

-0.09±0.20 

-0.17±0.23 

-0.30±0.23 

IR1, 46 

NGC 7027 

0.57±0.16 

0.57±0.18 

0.45±0.18 

0.22±0.19 

0.22±0.22 

0.10±0.22 

IR2, 47 


0.70±0.28 e 

0.70±0.30 e 

0.58±0.30 e 

0.18±0.22 

0.18±0.24 

0.06±0.24 

IR3, 47 


0.74±0.18 

0.74±0.21 

0.62±0.21 

0.19±0.23 

0.19±0.25 

0.07±0.25 

IR6, 47 

NGC 7354 

<0.87 

(< 0.93) 

< 0.45 

0.46±0.24 

(0.52±0.31) 

0.04±0.31 

IR1, 2 

NGC 7662 

<0.74 

< 0.87 

< 0.89 

0.32±0.21 e 

0.45d=0.23 e 

0.47±0.23 e 

IR2, 29 





0.14±0.68 

0.28±0.75 

0.30±0.75 

IR3, 29 

SwSt 1 

0.33±0.15 f 

0.58±0.29 





IR4, 11 

Vy 1-1 







IR1, 23 

Vy 1-2 

<0.85 

< 0.81 

< 0.78 

<-0.09 

< -0.13 

< -0.16 

IR1, 46 

Vy 2-2 

<-0.35 

< 0.33 d 

< -0.06 d 

-0.93±0.18 

-0.25±0.20 d 

-0.64±0.20 d 

IR1, 46 


-0.40±0.44 e 

0.28±0.44 d ’ e 

-0.11±0.44 d,e 

-0.89±0.31 e 

-0.21±0.31 d,e 

-0.60±0.31 d ’ 

e IR3, 46 


<-0.39 

< 0.29 d 

< -0.09 d 




IR4, 46 


a The [Kr/O] and [Se/O] values of Type I PNe are placed in parentheses, to indicate that O may not be a 
reliable metallicity indicator in these objects (see text). 

b Abun dances are relative to the lAsplund et al.l (120051) solar composition for ease of comparison with 
I Paper III . Using the lAsplund et alj d2Q09h solar composition, the listed [Kr/O] and [Se/O] values would 
increase by 0.05 and 0.02 dex, respectively, while [Kr/Ar] and [Se/Ar] would increase by 0.24 and 0.21 dex. 


c REFERENCES (O ptica l/UV/M 
(3) lAller et alj \ 1986); (4) lAller et 


Bernard-Salas et al. 

(I1999I) 


id-in frared Data ". . __ _ „ _ 

ahl dl988 ); (5) iBarken jl978allb|); (6) 


Alter fc Czvzald ( 19831); ( 2) lAller & Keves| |l987l) ; 


Bernard Salas et al 


(2003); (7) 


200^_(8) CpsuTetaLT 2004|/j_(9]_ ^uisi^Sr_et_alJ 1 19961); (10) IPe^TirccTfe Cro wtherl 

^ D^M^co_e^alJ_ 20fU h M2) lErcolan o et**aL _l 2004 ^_ (13) Feibelma^^^alJ jl99^) ; (14 ) 

iGirard e t al. I 2007^; (15) iGorny et all (1 20041) : (16 ) iHen rv et al]~i 200d. 12004 ); (17) 

(18) lHjrung_UAljer^ ^ (19lHvung etaL il993l ): (20) Hyungetal] jl994bl ); (21) 


Hyung et al. 


Hyung et al 


1995 ); 


1999a) ; 


fIR2nGebane"et al 


Kingsburgh & Barlow ( 

1994|); (27) Koppen et al. (1991); (28) Henrv et al. (i 200ol); Kwitter et al. 

200S) 

(29) Liu et al. (2004b|)a 

); (30) Lopez-Martin et al. (20021); (31) 

Henrv et alj ( 2000|^_ Milineo et al.l 

2002) 

(32) Perinotto & Corradi (199^); (33) Pottasch &: Beintema ( 

1999i); (34) Pottasch et al. t2000l) 

(35) 

Pottasch et al. 

(20051); (36) Pottasch & Surendiranath (200^); (37) Pottasch & Surendiranath ( 2007 

); (38) 
); (42) 

Pottasch et al. 

(2003); (39) Pottasch et al. (2003); (40) Ratag et alj il997f); (41) Samland et al.l (1992 

Shen et al. (20041); (43) 

Surendiranath et al. (2004|); (44) Tsamis et al. (20031); (45) Wesson &; Liu 

0 

0 

(M 


( 199ll); ( IR3) 


dl996h ; (IR6) iLumsdeiTet alj (12001 ); (IR7) llfaLd^^t alj d2001 ). 


Hora et al 


]l999l): ( IR4) 


(46) IWe sson et ahl (2005 ); (47) IZhang e t alj d2 005l). REFE RENCES (Nea r-Infrared Da ta); (IR1) Paperjij; 


Likkel et al.l d2006l) : (IR5) iLumsden &; Puxlevl 


d Ar is used as a reference element. This object exhibits N enrichments that fall just below the value for 
Type I classification, and hence its O abundance may not be reliable. 


e Based on a weak or uncertain detection. 


f ICF(Kr) from Equation [I] 

g IC 4997 has an uncertain O abundance, and Ar is used as the reference element for this object. 
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By compari ng our new analytical fits to the correlations plotted in Figured] (solid lines) 
with those from Paper I (dashed lines), it is expected that the ICFs and hence elemental Se 
and Kr abundances in PNe will be lower than previous estimates. Indeed, that is exactly 
what is seen (for comparison, see Table 10 of Paper II]). Overall, of the 33 PNe in which the 
Kr abundance could be determine d, t he average [Kr/(0,Ar)j = 0.82 ±0.29, compared to the 
mean of 0.98±0.31 found in Paper II . Likewise, the mea n [Se/(0,Ar)] = 0.1 2 ± 0.27 in the 


PNe exhibiting [Se IV] emission, down from Paper IPs 0.31±0.27. If the As plund et al. 


( 20091 ) solar abundances are used, the mean Se and Kr abundances in our sample would each 
increase by 0.05 dex. 

In individual objects, the derived Se and Kr abundances are typically lower by 0.1- 
0.3 dex compared to previous estimates, but_there are occasional exceptions. For example, 
the new Kr/Kr 2+ ICFs exceed those of Paper I when the S 2+ or Ar 2+ fractional abundances 
are sufficiently high. In the case of BD±30°3639, the Kr abundance actually inc reased b y 


0.14 dex from that found in iPaper III , since the S 2+ /S f raction o f 0.79 flBernard-Salas et al. 


20031 ) lies above the point where Equation Q] intersects Paper Js Equation 2 (see Figure [T] 
top left panel). For other PNe, the new Kr abundances do not exceed those from Paper II 
by more than 0.01 dex. In contrast, the Se abundance decreased in all objects of the sample. 
For some objects (e.g., Cn 3-1, M 1-11, and M 3-41) the Se abundances (or upper limits) 
are lower by as much as an order of magnitude. This is due to the extremely small 0 2+ /0 
fraction in these low-excitation objects, which led to a very large and uncertain ICF using 
Paper It s formula. For example, the Se ICF for M 1-11 decreased from 257±130 (IPaper III ) 
to 23.4±5.4 in the present study. In fact, the Se ICF formula from Paper I (their Equation 3) 
is negative and hence cannot be applied when 0 2+ /0 < 0.01626. Our Equation [8] is valid 
for all 0 2+ fractional abundances. This is the reason that upper limits to the Se abun dance 
are given in Table [17] for He 2-459 and M 1-12 when they were not determined in [Pa pe r II. 


The general reduction of Se and Kr abundances in this sample of PNe decreases the 
number of objects that are enriched by in situ s-process nucleosynthesis in their progenitor 
stars’ AGB stage of evolution. Since the sample is primarily composed of Galactic disk PNe, 
we follow [Paper III and assume that a PN is self-enriched in Se and Kr if its [Se/(O, Ar)] and 
[Kr/(0, Ar)] abundances exceed 0.3 dex. This cutoff lies above the dispersion of ~0.2 dex 
in the abundances of other light n-capture elements (Sr, Y, and Zr) in unevolved near¬ 
solar metallicity stars (IBurris et al.ll2000l: iTravaglio et ahll2004l ). and therefore (statistically) 
accounts for star-to-star scatter in initial Se and Kr abundances. Using this criterion, we 
find that 26 of 33 (79%) objects with determined Kr abundances are enriched by 5-process 
nucleosynthesis in their progenitor stars, while only 14 of 68 PNe (21%) are self-enriched in 
Se. Overall, 30 of the 81 PNe (37%) with determined Se and/or Kr abundances are 5-process 
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enriched^] The percentages of enriched PNe found in Paper II were 85% for Kr, 35% for Se, 
and 52% overall. 

In the 22 PNe in which both Se and Kr emission lines were detected , we find an abun¬ 
dance ratio [Kr/Se] = 0.5 ± 0.2 nearly identical to that derived in Paper II. This is in agree¬ 


ment (though toward the up per end of the range 


cleosynthesis in AG B stars (jClallino et al III 998 


with computationa 


Goriely & 


Mowlavi 


stuc 


2000 


ies_of s-process nu 


Busso et ah 2001 


Karakas et al. 2009 ). We refer the reader to Paper ij and Karakas et al. ( 2009 1 for a more 


extensive discussion comparing the observational results to predictions of nucleosynthesis 
models. 


5.1. Correlations Between Se and Kr Abundances and Other Nebular/Stellar 

Properties 


In [Paper ij . we searched for correlations between Se and Kr s-process enrichment factors 
and other nebular and stellar properties. We investigate here whether the improved Se and 
Kr _abundances affect those correlations. We do not repeat all details of the discussion from 
Paper II . and refer the reader to that paper for the full analysis and interpretation. 


One of the most significant correlations found in [Paper III was that PNe with higher 
mass progenitor stars (> 3-4 M 0 ; e.g., Type I and bipolar PNe) tend to show little if 
any s-pro cess enrichment of Se and Kr, compared to nebulae ejected by lower-mass stars. 
Foll owing iPaper III to impro ve statistics we consid e r a P N to be Type I if it meets either 
the iPeimbertl (119781) or the Kingsburgh fe Barlowi (119941 1 criteria. Table [TH] displays the 


average [Se/(0, Ar)] and [Kr/(0, Ar)] abundances for various subclasses of PNe, including 
Peimbert types, different morphologies, and central star types. The average Se abundance is 
0.39 dex higher in non-Type I PNe than in Type I objects, while the average Kr abundance 
is nearly an order of magnitude higher. Indeed, no Type I PN in this sample is enriched in 


eithe r Se or Kr by more than 0.05 dex (0.3 de x using the solar abundances of lAsplund et al. 


2009). Kolmogorov-Smirnov (K-S) tests (e.g., Press et al. 1992) indicate that the differences 


in s-process enrichments in Type I and non-Type I PNe are statistically significant, with 
probabilities p^ s = 0.025 and 0.009 that the Se and Kr abundances (respectively) in these 
two subclasses are drawn from the same distribution function. This result is robust against 
the metallicity reference element (Ar or O) used for non-Type I PNe, as well as the adopted 
solar composition. 


2 The percentages of enriched objects increase by less than 5% when adopting the Asolund et al. ( 2009 1 


solar abundances. 
























































Table 18. Kr and Se Abundances Vs. Nebular and Stellar Properties 



Mean 

Mean 


Number of 

Mean 

Mean 


Number of 

Property 

[Se/(0, Ar)] a 

[Se/(O, Ar)] b 

< <7 > C 

Se Detections 

[Kr/(O, Ar)] a 

[Kr/(O, Ar)] b 

< a > c 

Kr Detections 

Type I 

-0.22 

-0.01 

0.27 

12 

-0.11 

0.13 

0.12 

3 

Non-Type I 

0.17 

0.20 

0.25 

55 

0.85 

0.91 

0.25 

30 

Bipolar 

0.08 

0.13 

0.38 

14 

0.50 

0.58 

0.28 

8 

Elliptical 

0.15 

0.19 

0.24 

28 

0.95 

1.00 

0.36 

15 

[WC] 

0.19 

0.24 

0.25 

16 

0.79 

0.84 

0.31 

10 

WELS 

0.09 

0.12 

0.24 

12 

0.45 

0.50 

0.16 

3 

Non-[WC]/WELS 

0.10 

0.15 

0.28 

39 

0.86 

0.92 

0.30 

20 

Full Sample 

0.12 

0.17 

0.27 

67 

0.82 

0.87 

0.29 

33 


Note. — Only PNe exhibiting Se and/or Kr emission and with determined O and Ar abundances are included. 

a Ar is used as the reference element for Type I PNe, as well as M 1-31, NGC 688 1, and Vy 2-2 (w hich have N/O ratios just below 
the Type I cutoff); O is used for all other objects. Solar abundances are taken from Asplund et al. (120051 ). 

b Solar abundances from Asplund et al. ( 2009f) . 


CT 

GO 


The < er > are the mean absolute deviations in the Se and Kr abundances. 
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Bipolar PNe are believed to derive from a more massive population than other mor- 


--—-- — — —— 

and tendency to exhibit Type I compositions ( 

Cornu 

i & Schwarz 

1995; Gornv et al. 

1997; 

Torres-Pcimbert & Peimbert 

1997; 

Stanghcllini et al. 

2002). The average [Se/(0, Ar)] and 


[Kr/(0, Ar)] abundances in bipolar PNe are lower than in elliptical PNe (Table [18]), though 
the difference is not as striking as for Type 1 and non-Type 1 objects. Unlike Type 1 PNe, 
several bipolar objects show evidence for self-enrichment in Se and Kr. K-S tests suggest 
that the Se and Kr abundances in bipolar and elliptical PNe are not strongly distinct, with 
respective probabilities 0.42 and 0.21 that the abundances are drawn from the same distribu¬ 
tion function. Paper II found a similar result, suggesting that morphology is not as reliable 
as chemical composition in gauging the masses of PN progenitor stars. 


We also confirm Paper Ill's result that there is no significant difference in the s-process 
enrichments of PNe with H-rich and H-poor ([WC] and weak emission line stars, or WELS) 
central stars. [WC] and non-[WC] PNe have similar average Se and Kr abundances, and K-S 
tests indicate that the Se and Kr abundance distributions are the same in the two subclasses 
(Pfcs(Se) = 0.79 and pfc s (Kr) = 0.80). The average enrichment factors are smaller for PNe 
with WELS nuclei, but K-S tests show that there is little evidence that the Se and Kr 
abundance distributions are distinct from those of PNe with other central star types. These 
results adhere to previous findings that the nebular compositions of [ WC] and WELS PNe 


do not show significant differences from those with H-rich nuc l ei flGornv &; Stasinskal 11995 


Pena et al.l l200ll; iDe Marco fe Barlowll200ll; iGirard et al.l l2007t iGornv et al.l 120091) . 


Finally, in iPaper II we found that Se and Kr enrichments positively correlate with 


nebular C/O ratios, as is e xpected theoretically (e.g.. iGallino et al.lll998l: iBusso et al.l 12001 


Karakas fe Lattanzio 


20 


stars flSmith fe Lambert 


4f) and seen for other n-capture element s in AGB and post-AGB 


1990l; lAbia et al.ll2002l; Ivan Winckell 120031 ). In Figure El we plot the 


new Se and Kr abundances against log(C/ Q), where the carbon abundances were determined 
from UV forbidden lines. As discussed in Paper II, Hb 12 is an outlier in the [Kr/(0, Ar)] 
plot in that it exhibits a high Kr abundance but a low C/O ratio. We have excluded Hb 12 
from this analysis, but note that its inclusion greatly weakens the correlation between Kr 
and C abundances. 

The correlation coefficients r are given in each panel of Figure El along with the prob¬ 
ability of no correlation p r= o (i.e., the probability that the observed correlation between the 
same number of objects could have been found from two uncorrelated populations). The 
values of r indicate marginal correlations (slightly less robust than found in Pap er II, par¬ 
ticularly for Kr versus C/O), likely due to the relatively large observational uncertainties in 
both the Se and Kr abundances and in the C/O ratios. 
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Fig. 6.— [Se/(0, Ar)] (left) and [Kr/(0, Ar)] (right) are plotted against C/O (derived from 
UV forbidden lines) for PNe in the sample of lPaper III . Typical error bars are shown in the 
left panel. The solid lines indicate the best linear fits to the correlations. The correlation 
coefficient r and significance p v=0 are given in each panel. The discrepant point in the right 
panel corresponds to Hb 12. 


We applied a least-squares fitting routine to find the best linear fit to each correlation. 
These are plotted as solid lines in Figure [6], and correspond to the equations: 


and 


[Se/(0, Ar)] = (-0.01 ± 0.05) + (0.40 ± 0.16)log(C/O), 


[Kr/(0, Ar)] = (0.30 ± 0.11) + (0.69 ± 0.32)log(C/O), 


( 10 ) 


( 11 ) 


where the uncertainties are derived from the dispersion in the fits. The fits have lower inter¬ 
cept values than Equations 4 and 5 of Paper II, but the slopes are quite similar. Therefore, 
while the correlations between Se and Kr abundances and the C/O ratio of PNe are slightly 
less robust than found in lPaper II . the abundances follow nearly the same trend. 


We have therefore verified the correlations found in [Paper III , in particular that Type I 
PNe show little to no s-process enrichment compared to their counterparts with less massive 
progenitor stars, and that there is evidence for a positive correlation between Se and Kr 
enrichments and nebular C/O ratios. In addition, we confirm that there is no evidence that 
Se and Kr enrichment factors differ in PNe with H-rich nuclei and those with [WC] or WELS 
central stars. The results do not depend on whether Ar or O is used as the reference element 
in non-Type I PNe, nor on the adopted solar composition. Thus, the new ICFs for Se and 
Kr have corrected a systematic overestimation of Se and Kr abundances by 0.1-0.3 dex for 
most PNe, but the newly-derived abundances do not alter the overarching conclusions of 
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Paper II. 


5.2. 5-Process Enhancements for Thick Versus Thin Disk Nebulae 


In order to assess whether a PN is self-enriched in Se or Kr, and by what amount, 
we must adopt initial abundances for these species in the progenitor star and compare 
them to the observed values. Since the initial abundances cannot be determined directly, 
we estimate them from abundances of measurable species and the characteristics of the 
Galactic population to which the progenitor star belongs. Thus far we have used 0 and Ar 
as proxies for the overall metallicities of PN progenitor stars, and taken the initial values of 
Se/(0, Ar) and Kr/(0, Ar) to be solar. This assumption, implicit in our previous criterion 
for designating a PN as self-enriched, is not necessarily valid if the sample includes objects 
from multiple stellar populations. 


In this section, we relax the assumption that the PN progenitor stars all have scaled 
solar compositions, and take into account the well-established abundance patterns in the 
subpopulations of the Milky Way disk. Many large spectroscopic surveys have attempted to 
characterize and distinguish the thin and thick disks in t erms of properties such as kinemat- 


i . 

Bensbv et al. 

200j 

2004; 

Reddv et a' 

. 2006 

Feltzing & Chibk 

2013; Bensby et al. 2 

03). 


Ramfrez et al 


] 2007t IaI 


ves-Brito et ah 


2002 


2010 


metallicity ([Fe/H]) distribution functions (MDFs) that overlap heavily, but all studies find 
that the thick disk’s MDF is skewed toward lower [Fe/H] values than the thin disk. The 
centroids of the M DFs are not always wel l -defined and may vary with the characteristics 
of the sample (e.g., ICasagrande et al.l 1201 ll: iKordopatis et al.l l201lh . but are approximately 
< [Fe/H] >~ —0.6 ± 0.1 dex for the thick disk and < [Fe/H] —0.1 ± 0.1 dex for the thin 
disk. 


In view of the large spread in [Fe/H] within each disk component, it would be inappro¬ 
priate to assume a single characteristic value for PNe of each population. However, a robust 
and much tighter relation emerging from these spectroscopic surveys is the separation of the 
two populations in their [a/Fe] ratios, with thick dis k values M).2 0.8 dex higher than the 


thin disk for moderately subsolar metallicities (e.g., iReddv fe Lambertl 120091; lAnders et al. 


2014; Bensby et ah 2014). While there is controversy over whether this quantity is essentially 


a step function with a break between the two populations or two parallel linear correlations 
separated vertically in an [a/Fe] vs. [Fe/H] diagram, we will use the former approximation 
below for simplicity. (Note that this phenomenon introduces an ambiguity in the use of the 
a elements O and Ar as proxies for [Fe/H] in PNe: if the sample includes both thin and 
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thick disk objects, [Fe/H] cannot be unambiguously inferred from [ct/Fe] because [a/Fe] is 
double-valued.) 


In Paper II, we examined correlations between Se and Kr abundances and several char¬ 
acteristics of the PNe in onr sample, but did not attempt to distinguish between different 
Peimbertl (119781 ) classes, except to pay special attention to objects of Type I, which de¬ 
scend from relatively high-mass progenitors. Peimbert also defined two additional classes of 
Galactic disk PNe, Types II and III. Type III PNe have radial velocities |u ra d| > 60 km s -1 , 
and in terms of kinematics represent an older stellar population than Type II PNe, which 
exhibit lower radial velocities. In their spatio-kinematic characteristics, Type II PNe are 
now r ecogni z ed as members of the thin disk, while Type III PNe belong to the thick disk 


Phillips 2005 )1^1 The O, Ne, S, and Ar abunda nces o f Type III PNe are only slightly 


(e.g., 

lower than in their Type II counterparts (IPerinotto et al.l I2004J: IStanghellini fc Havwood 
20101 ). However, since these elements are a species, this does not preclude the possibil¬ 
i ty that Type III PNe are deficient in Fe-pea k elements. Measurements of [Zn/H] in PNe 
(IDinerstein fe Geballell200ll; ISmith et ahll2014l Dinerstein et al. 2015, in preparation), an ef¬ 
fective proxy for [Fe/H], demonstrate that many high-velocity (Type III) PNe have subsolar 
metallicities. The a-enhancement at lower [Fe/H] values thus partially compensates for low 
Fe abundances, yielding similar a-element abundances for Type II and Type III PNe. 

We separate the full sample of Paper II into two subsets, Type II and Type III PNe@ 
which at least statistically correspond to the thin and thick disks, respectively. To avoid 
circular reasoning, we base our population assignments purely on kinematical grounds rather 
than abundances. This criterion is not ideal, since kinematic mixing blurs the boundaries 
between the thin an d thick disk popula t ions even when all t hree kinematic components (U, 
V, W) are known (lAnders et al.l 12014 iBensbv et al.l 120141) . but for PNe radial velocities 
are generally the only_available kinematic information. We use the radial velocities of the 
Strasbourg catalog (lAcker et al.l 1199411 . which provides data for 109 objects in our sample, 
and apply Peimbert’s original criterion of |u ra d| > 60 km s -1 to distinguish Type III (thick 
disk) from Type II (thin disk) PNe. Given large uncertainties due to the absence of proper 
motion data and poorly-known distances, we do not attempt to correct the heliocentric 
velocities for solar motion or Galactic rotation. (Type III and Type II ob jects app ear to 
be kinematically distinguishable without this correction; see Fig. 2 of Phillips 20051) . This 


3 While the modern nomenclature of thick and thin disks had not yet been developed at the time, P eimb ert 


((l978h proposed that in view of the similar kinematics of Type III PNe to those of old metal-deficient stars, 
the PNe probably were also Fe-poor, although no means was available to verify this suggestion. 

4 Both groups contain several nominally N-rich Type I PNe. The inclusion of these objects reduces the 
average abundances discussed below by 0.03-0.05 dex, but does not affect our conclusions. 









































yields 39 PNe in the high-velocity Type III group, and 70 in the Type II group, including 
objects for which only upper limits to the Se and/or Kr abundances are available. 


We extrapolate each PN progenitor star’s elemental Fe abundance using its measured 
[Ar/H] abundance and the average [q/Fe] of the disk population to which the object is 
kinematically associated. We use Ar as the a tracer for all sources, to avoid possible system¬ 
atic effects; stellar abundance studies tend to find a larger enhancement for [O/Fe] than for 
heavier species such as Si, Ca, and Ti (e.g., Bensbv et al. 2014). While this may be partly 
due to systematic uncerta inties caused by using different spectros copic indicators to deter¬ 
mine O abundances (e.g., iBensbv et al.ll2004r iRamfrez et al.l 120071) . the greater enrichment 
of lighter q-elements is also likely a consequence of metallicity-dependent yields in massive 
stars (e.g., iMcWillia/m et ah 2008). Our choice of [Ar/H] also avoids possible artificial effects 
introduced by using different q tracers for different PNe. Although there are no direct data 
on the abundance of Ar in stars, we take typical average values for the heavier a species in 
the metallicity regim e represented in ou r sample: [q/Fe] = 0.25 dex for the thick disk, and 
0.0 for the thin disk ( Anders et al. 2014 : Bensbv et ah 2014). 


The next step is to extrapolate, from the inferred Fe abundances, the initial abundances 
of Se and Kr in the PN progenitor stars. These two elements fall in a comparative desert 
of stellar abundance data in the periodic table, between the high-mass tail of the iron peak 
at Zn and the widely-observed light -s peak of Sr, Y, and Zr. A few direct measurements of 


Se abu ndances in stars have recently been obtained by iRoedererl (120121) and iRoederer et ah 


(2014). These authors found a roughly constant value of [Se/Fe] = 0.16 ± 0.09 for seven 
stars spanning —2.8 < [Fe/H] < —0.6. We adopt this value for our Type III (thick disk) 
subgroup, and assume a solar Se/Fe ratio for the Type II (thin disk) objects. By applying 
the [q/Fe] and [Se/Fe] values for the two groups to the measured [Ar/H] abundances, we 
infer the initial [Se/H], The net effect is [Se/H] = [Ar/H] — 0.09 dex for the thick disk and 
[Se/H] = [Ar/H] for the thin disk. 

Since Kr cannot be observed in cool stars, we approximate its behavior relative to Fe 
by considering the ele ments bracketing it in the periodic table. In addition to their results 
for Se, Roederer et al.J (2014) found a slightly higher average abundance for As than for Se: 
[As/Fe] = 0.28 ± 0.14, and no trend with [Fe/H]. [Sr/Fe] is slightly elevated (0.1-0.2 dex) 
in unevolved metal-poor stars (llshigaki et al.l 120131 ). wh ile Y abundances are fa irly constant 
and slightly subsolar in disk stars with [Fe/H] > —1 (iMishenina et all 201 31). Given this 
information, we believe that it is reasonable to adopt the same initial values for [Kr/Fe] (and 
hence [Kr/H]) as for [Se/Fe], 


The enrichment factors for each nebula are computed from the difference between the 
measured values of [(Se, Kr)/H] and the initial values predicted by the procedure described 
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above. Se was detected in 42 of the 70 low-velocity PNe and in 23 of the 39 high-velocity 
objects. We find that the average self-enrichment of Se is 0.07±0.39 dex in thin disk PNe and 
0.21±0.25 dex in thick disk PNe. In contrast, the average enrichment of Kr is considerably 
higher in both populations: 0.62±0.41 dex in 23 thin disk PNe and 0.75±0.33 dex in seven 
thick disk objects. Note that the numbers after the ± symbols result at least in part from 
intrinsic dispersion in the actual range of enrichments, and should not be interpreted simply 
as errors in the average values. 

These results suggest that thick disk PNe may exhibit larger (by 0.1-0.2 dex) Se and 
Kr enrichment factors on average than their thin disk counterparts. Our assumption that 
the thick disk members have lower initial Se and Kr abundances than thin disk PNe may 
contribute to this result, but does not fully account for it. Moreover, there is strong evidence 
for systematic differences in abundances between the thick and thin disks from stellar surveys, 
and this is supported by the significantly subsolar Se abundances in some of our observed 
PNe (presumably from stars that experienced little or no s-processing). 


The possible offsets in Se and Kr self-enrichments are in qualitative agreement with the¬ 
oretical calculations of s-process yields for AG B stars that predi c t a metallicity depen dence 
for the production of the light-s species (e.g., Busso et al.l 12001 : ICrist alio et al.l 120091) . We 


examined this effect in the FRUITY database j fCristallo et al. 


2011). For stars of initial 


mass 1.5-2.0 M 0 the Se s-process enhancement increases by 0.04-0.08 dex when going from 
[Fe/H] — 0 (Z — 0.014) to [Fe/H] = —0.67, and by 0.24 dex at [Fe/H] = —1.15. For Kr 
the corresponding values are 0.09-0.15 dex larger at [Fe/H] = —0.67 and 0.3 dex larger at 
[Fe/H] = —1.15 compared to solar metallicity models. These models also predict higher Se 
and Kr yields for 2.5-3.0 M 0 stars as metallicity decreases. However, higher-mass stars may 
have initially higher rotational velocities, which results in lower s-process production rates 
flPiersanti et al.l 20131 ). At lower metallicities ([Fe/H] < —1), there is large scatter in the 
observed (initial) abundances of light-s species in unevolv ed stars, but this is mainly due 
to the operation of mechanisms other than the s-process (jSneden et ah 2008). Our sample 
is drawn from more metal-rich populations and is therefore unlikely to include a significant 
number of such low-metallicity objects. 


While we have found tentative evidence for a metallicity dependence in the enrichments 
of Se and Kr, these enhancements also depend on initial stellar mass. Additional factors 
such as rotation and details of convection may also play a role. Differences in average values 
between different populations will therefore depend on the distribution of stellar masses as 
well as metallicities, leading to a natural dispersion in the enrichment factors within each 
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population. This expectation is supported by the large standard deviations in the average 
Se and Kr enhancements in thin and thick disk PNe, which indicate a substantial spread 
and overlap in the populations. It would be interesting to extend this kind of comparison to 
additional s-process products, for example nuclei belonging to the heavy-5 peak. It will also 
be useful to develop more robust methods for determining the parent population of individual 
PNe, such as obtai ning additional kinematic in f ormation and/or F e-group abundances using 
[Zn/H] as a tracer ( Dinerstein fe Geballe 2001 : Smith et al. 2014 1. 


6. CONCLUSIONS 


We have conducted a study of the ionization balance of the n-capture elements Se 
and Kr in ionized nebulae, utilizing recently-determined atomic data for photoionization, 
recombination, and charge exchange processes for relevant ions of these elements. The new 
atomic parameters were incorporated into the C13 release of Cloudy OFerland e t al. 20131), 
which we used to derive new Se and Kr ionization correction factor formulae. We app lied 
these ICF prescriptions to our near-infrared survey of [Kr III] and [Se IV] lines (Paper 11) to 
derive revised Se and Kr abundances in a sample of 120 PNe. Our main conclusions are: 


1. We tested the ability of models that incorporate the new atomic data to reproduce 
the observed Kr ionization balance in 15 PNe that exhibit emission lines from multiple 
Kr ions. We found a systematic disagreement between the modeled and observed Kr 
ionization balance for all 15 PNe, regardless of stellar and nebular properties. The 
[Kr III] intensities were predicted by the models to be stronger than observed, while 
the [Kr IV] lines were predicted to be too weak. Given that no systematic discrepancies 
were seen in the ionization balance of other elements, the most probable explanation 
is that the Kr photoionization and recombination data are inaccurate. To correct 
this disagreement, we adjusted the photoionization cross sections of Kr + -Kr 3+ within 
their stated uncertainties, and altered the DR rate coefficients for the same ions by an 
amount somewhat larger than the approximate uncertainties cited by Ste rling (20111). 
Because only one Se ion has been unambiguously detected in astrophysical nebulae, a 
similar test of the Se atomic data is not currently possible. 


This result underlines the importance of the interplay between theoretical and ex¬ 
perimental atomic physics and observational spectroscopy. All three approaches were 
needed to best reproduce the observed Kr ionization balance in PNe. This provides 
added motivation for deep, high-resolution optical and infrared spectroscopy of PNe 
to detect multiple ionization states of n-capture elements. Such observations enable 
more accurate n-capture element abundances to be derived and provide critical tests 
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to newly-determined atomic data, both of which are needed to advance the efficacy of 
nebular spectroscopy for the study of trans-iron element nucleosynthesis. Our study 
also highlights the need for more accurate (e.g., experimental) DR measurements for 
low-charge ions of heavy elements. In a forthcoming paper, we will explore the sensi¬ 
tivity of Se and Kr abundance determinations to uncertainties in photoionization and 
recombination data. This will help to elucidate which ionic systems are most in need 
of further investigation. 


2. We computed multiple grids of Cloudy models, spanning a wide range of central star 
temperatures, luminosities, nebular densities, and chemical compositions in order to 
derive new ionization correction schemes for unobserved Se and Kr ions. ICF formu¬ 
lae were computed from correlations between Se and Kr ionic fractions and those of 
commonly-detected ions of lighter elements. We determined analytical ICFs for single 
Se and Kr ions and for combinations of those species, which can be applied to optical 
and near-infrared observations of ionized nebulae. We tested the various ICF formulae 
with the optical and near-infrared spectra of the 15 individually modeled PNe. The 
ICF-derived Se and Kr abundances agree within the uncertainties with the modeled 
values in nearly all cases. The abundances determined with ICF formulae that rely on 
single Kr ions show larger scatter than those computed from multiple ionization states, 
as is expected. 


3. The new ICFs were applied to the near-infrared survey of [Paper III to re-evaluate Kr 
and Se abundances in a sample of 120 PNe. For most objects, the revised Se and 
Kr abundances are 0.1-0.3 dex lower than previous estimates. We derive new average 
abundances of [Se/(0, Ar)] = 0.12 ± 0.27 and [Kr/(0, Ar)] = 0.82 ± 0.29 in PNe in 
which Se and Kr emission lines were detected , compared to the respective values of 
0.31 ± 0.27 and 0.98 ± 0.31 found in Paper IF These average abundances show that 
Se is not enriched by in situ s-process nucleosynthesis in many of the PNe in which 
it was detected. Indeed, according to our criterion that requires [(Se)/(0, Ar)] and 
[Kr/(0, Ar)] > 0.3 dex to be considered self-enriched, only 21% of the PNe exhibiting 
Se emission are enhanced in this element. Many PNe have subsolar Se abundances, 
which may reflect their initial compositions. In contrast, Kr abundances are enhanced 
relative to solar in 79% of the PNe in which [Kr III] was detected. 


4. While the overall Se and Kr abundances in the sample of Paper II decreased from 
their previous values, we confirm the correlations found in that paper between Se and 
Kr abundances and other nebular and stellar properties are unaffected. Specifically, 
Type I PNe show no evidence of significant s-process enrichments, and bipolar nebulae 
tend to have lower Se and Kr abundances than elliptical PNe. Since Type I and 
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bipolar PNe tend to be associated with more massive progenitor stars (> 3-4 M 0 ), this 
suggests that intermediate-mass AGB stars experience smaller s-process enrichments 
than lower-mass stars, at least in the case of Se and Kr. In addition, the positive 
correlation between Se and Kr abundances and the C/O ratio found in Paper II is seen 
with the newly-derived abundances. Finally, we do not find significant differences in 
the s-process enrichment of PNe with different central star types ([WC], WELS, and 
H-rich). These results do not depend on the choice of the metallicity reference element 
for non-Type I PNe or on the adopted solar composition. 

5. We took into account the initial compositions of the PNe by dividing the sample into 
two subgroups based on their radial velocities: low-velocity (Peimbert Type II) and 
high-velocity (|n ra d| > 60 km s -1 ; Peimbert Type III). We argue that these subgroups 
statistically correspond to the Galactic thin and thick disk stellar populations, respec¬ 
tively. To estimate the initial Se and Kr abundances of PNe in each subgroup, we 
utilize the [Se/Fe] abundances found in metal-poor and near-solar metallicity stars. 
The initial Se and Kr abundances were inferred using the abundance of the a el¬ 
ement Ar to indicate the initial composition of each nebula, taking into account 
the differing [a/Fe] abundances in thin and thick disk populations. We find that 
< [Se/H] >= 0.07±0.39 dex in the thin disk sample and 0.21±0.25 dex in thick disk 
objects. A similar trend is seen for the average Kr enrichments in thin and thick disk 
PNe: < [Kr/H] >= 0.62±0.41 and 0.75±0.33 dex, re spec tively. Therefore, when ac¬ 
counting for the initial compositions of the PNe in the Paper IJ] sample, the average Se 
and Kr enrichment factors appear to be larger by 0.1-0.2 dex in metal-poor thick disk 
objects than in thin disk PNe^ This o ffset is in qualitative agreement with models of 
AGB nucleosynthesis (e.g., Cristallp et al, 2011) for 1.5-2.0 M 0 stars with metallicities 


— 1.15 < [Fe/H] < 0. To further investigate the possible offset in s-process enrichments 
in thin and thick disk PNe, another method is needed to validate the association of the 
high-velocity PNe with the thick disk, for example Zn abundance determinations (a 
proxy for Fe abundances). It would also be of interest to study the behavior of other 
n-capture elements in PNe from the different disk populations. 
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